o 

(N 






Understanding the J/ip Production Mechanism at PHENIX 

by 

Todd Kempel 

A dissertation submitted to the graduate faculty 

in partial fulfillment of the requirements for the degree of 

DOCTOR OF PHILOSOPHY 



^ Major: Nuclear Physics 

I 

O 

^ Program of Study Committee: 

Ph John G. Lajoie, Major Professor 

Kevin L De Laplante 
^"^ Soren A. Prell 

J^ Jorg Schmalian 

Q\ Kirill Tuchin 

(N 

o 



Iowa State University 

Ames, Iowa 

2011 

Copyright @ Todd Kempel, 2011. All rights reserved. 



TABLE OF CONTENTS 



LIST OF TABLES m 

LIST OF FIGURES EEl 

CHAPTER 1. Overview [D 

CHAPTER 2. Quantum Chromodynamics [3] 

2.1 The Standard Model [3] 

2.2 Quarks and Gluons [5] 

2.3 Asymptotic Freedom and Confinement [6] 

CHAPTER 3. The Proton E 

3.1 Cross-Sections and Luminosities [5] 

3.2 Deep-Inelastic Scattering [10] 

3.3 Structure Functions and Bjorken Scaling [12] 

3.4 Altarelli-Parisi Evolution [13] 

CHAPTER 4. Transversely Polarized Protons [IT] 

4.1 Experimental Observations [18] 

4.2 The Sivers Effect [21] 

4.3 Transversity [23] 

4.4 The Boer-Mulders Effect [21 

4.5 The Coffins Effect M 

4.6 PDF Extractions and Global Fits [23 

CHAPTER 5. The J/ip Meson M 

5.1 Charmonium [31] 



HI 

5.2 The OZI Rule [32] 

5.3 Angular Distributions of Decay Products [33] 

5.4 Reference Frames [37] 

5.5 Production Mechanism [5T] 

5.5.1 Color Evaporation Model [U 

5.5.2 Color-Singlet Model \M 

5.5.3 Color-Octet Model M 

5.6 Discussion [l5] 

CHAPTER 6. RHIC and the PHENIX Experiment M 

6.1 RHIC M 

6.2 PHENIX M 

6.2.1 Detectors for Event Characterization [52] 

6.2.2 Local Polarimeter [52] 

6.2.3 Central Arm Detectors [53] 

6.2.4 Muon Arm Detectors \^ 

6.2.5 Triggering [58] 

6.2.6 J/V' mesons at PHENIX [60] 

CHAPTER 7. J/V" Transverse Single Spin Asymmetry [62] 

7.1 Motivation [62] 

7.2 J/V' Selection [63] 

7.3 Asymmetry Formulae [67] 

7.4 Polarization and Geometric Scale Factors [73] 

7.5 Systematic Uncertainties [71] 

7.6 Results and Discussion [75] 

CHAPTER 8. J/ip Angular Decay Coefficients [80] 

8.1 Introduction [SU] 

8.1.1 Measuring A^, A^, and A^^ [80] 

8.1.2 Simulation Procedure [82] 



IV 

8.2 Quality Assurance [83] 

8.2.1 Run Selection M 

8.2.2 J/V' Selection [85] 

8.3 Low Level Simulation Tuning [86] 

8.4 J/V' Yield M 

8.5 J/V' Angular Decay Coefficients [97] 

8.6 Discussion [98] 

CHAPTER 9. Conclusions and Future Measurements [TU51 

BIBLIOGRAPHY [TU^ 

APPENDIX A. Background Contributions to Angular Distributions .... 11191 

APPENDIX B. Relative Acceptance for Decay Angles \TM 

APPENDIX C. Data Tables [I22] 



LIST OF TABLES 



Table 2.1 Properties of the Standard Model gauge Bosons U] 

Table 2.2 Properties of the Standard Model quarks and leptons H] 

Table 5.1 Dominant J /il) decay modes [3H 

Table 6.1 Energy, spatial, and timing resolution of the PHENIX PbGl and PbSc 

calorimeters for electrons [5^ 

Table 7.1 Background fractions for J /ip mesons in the central and muon spec- 
trometers at several transverse momenta [68l 

Table 7.2 Background asymmetries as a function of pT for the PHENIX muon 

spectrometers [72l 

Table 7.3 Geometric scale factors for the central spectrometer [7^ 

Table 8.1 Frame-invariant angular decay coefficient. A, calculated for a single in- 
tegrated data point 11021 

Table C.l J/ip SSA vs. pT in forward, backward and central rapidities 11231 

Table C.2 J /ij: meson angular decay coefficients and uncertainties binned in pT . 11241 
Table C.3 Average x^ /ndf for two-dimensional fits in each arm to determine an- 
gular decay coefficients as a function of pr- Each fit has 36 degrees of 

freedom 11251 

Table C.4 J /ij) meson angular decay coefficients and uncertainties binned in xp ■ 11261 



VI 



Table C.5 Average x^ /ndf for two-dimensional fits in each arm to determine an- 
gular decay coefficients as a function of xp- Each fit has 36 degrees of 
freedom [HT] 



vu 



LIST OF FIGURES 



Figure 3.1 Scattering of a 188 MeV electron beam on a hydrogen target [TT] 

Figure 3.2 Commonly used variables in DIS [T^ 

Figure 3.3 Compilation of measurements of the F2 structure functions [T5l 

Figure 3.4 Parton Distribution Functions from the MSTW collaboration calculated 

at next-to-leading order in a^ [T6l 

Figure 4.1 Fermilab E704 transverse SSAs for tt^, -k^ and vr^ from p+p collisions 

at -y/s = 20 GeV as a function of Feynman-x (from (5)) [20l 

Figure 4.2 Brahms transverse SSAs for vr^ and vr^ from p+p collisions at ^/s = 
62.4 GeV as a function of Feynman-x for various ranges of transverse 

momenta (from (42)) [^ 

Figure 4.3 STAR transverse SSAs for vr'^ from p+p collisions at y^ = 200 GeV as 

a function of Feynman-x (from (2)) [2T] 

Figure 4.4 Sivers function from a global fit to SIDIS data [26l 

Figure 4.5 Preliminary PHENIX transverse SSAs for vr" from p+p collisions at 

^/s = 200 GeV as a function oi pT [271 

Figure 4.6 Incoming partons from processes contributing to vr" production at mid- 
rapidity from p+p collisions at y^=200 GeV US] 

Figure 4.7 Transversity distribution from a global fit to SIDIS and e^e~ data. . . [29l 

Figure 4.8 Collins function from a global fit to SIDIS and e'^e" data [30l 

Figure 4.9 Boer-Mulders function from a proof-of-principle global fit to SIDIS data 
assuming that the Boer-Mulders function is proportional to the Sivers 
function [30] 



Vlll 



Figure 5.1 Various charmonium resonances and their decay channels [321 

Figure 5.2 Charmonium resonances confirmed in data from various experiments. . [33l 

Figure 5.3 Example of a decay process suppressed by the OZI rule [331 

Figure 5.4 Angular distributions of decay products from a spin-1 boson decaying 

into two spin-1/2 fermions [551 

Figure 5.5 Illustration of i-axes, i?, and (/? in several reference frames I39l 

Figure 5.6 Orientation of the y-axis depending on the direction of the boost. . . . I39l 
Figure 5.7 Comparison of the Color Evaporation Model at NLO in as with J/ip 

cross-sections from CDF H3l 

Figure 5.8 Comparison of the CSM at NLO in as with J ftp cross-section data from 

RHIC M 

Figure 5.9 Comparison of the CSM at NLO in as with J/ip spin- alignment data 

from PHENIX M 

Figure 5.10 Feynman Diagrams for color-singlet Charmonium Production B51 

Figure 5.11 Feynman Diagrams for color-singlet Charmonium Production H6l 

Figure 5.12 Comparison of the COM at NLO in as with the J/ip cross-section from 

CDF M 

Figure 5.13 Comparison of the COM at NLO in as with J/ip spin- alignment data 

from CDF M 

Figure 6.1 Setup of the RHIC hydrogen-jet polarimeter [^ 

Figure 6.2 Setup of the RHIC carbon-target polarimeter [STl 

Figure 6.3 PHENIX Detectors in beam and side views |5T1 

Figure 6.4 Diagrammatic view of a PHENIX Pb-scintillator calorimeter module. . I55l 

Figure 6.5 Exploded view of a PHENIX Pb- glass calorimeter super- module. . . . I56l 

Figure 6.6 PHENIX Pad Chamber pad and cell configurations 1571 

Figure 6.7 Cut-away view of the PHENIX muon detectors [551 

Figure 6.8 PHENIX magnetic field orientation during 2009 RHIC run [59l 

Figure 6.9 Logic for a 'deep' ;uld LLl trigger [60] 



IX 



Figure 6.10 Logic for a 'shallow' //Id LLl trigger [6T] 

Figure 7.1 Transverse single-spin asymmetry generating diagrams for J/ip meson 

production [63l 

Figure 7.2 Invariant mass spectrum for electron candidate pairs (solid circles) and 

uncorrelated track pairs (open circles) in the central spectrometer. . . [65l 

Figure 7.3 Invariant mass spectrum for oppositely-charged muon candidate pairs 

and charged candidate pairs with the same sign [66] 

Figure 7.4 Invariant mass spectrum of electron pairs showing the contributions 

from all physical sources [671 

Figure 7.5 Number of e^e~ pairs with invariant masses Me [2.7, 3.4] GeV/c^ di- 
vided by number of counts in the BBC integrated over individual stores, 
plotted against store number, and fit with a constant [69l 

Figure 7.6 Kinematics of accepted e+e" pairs for J/ijj An analysis [701 

Figure 7.7 J/^ transverse SSA measured at PHENIX from ^/s =200GeV p+p 

collisions plotted against transverse momentum [771 

Figure 7.8 J /ij: transverse SSA measured at PHENIX from ^/s =200GeV p-\-p 

collisions plotted against Feynman-x [78l 

Figure 7.9 x range covered for various pT ranges of vr'^ production at mid-rapidity in 
the PHENIX central spectrometers and for the parton in the polarized 
proton from J /ip production in the PHENIX muon spectrometers at 
xf>^ from p+p collisions at -v/s=200 GeV [79l 

Figure 8.1 pz at the first /iTr station from a full simulation with A^=l in the 
Helicity frame divided by a full simulation with A^=0 for the North 
(left) and South (right). The A^=l distribution is clearly shifted to 
larger p^ Dashed lines show the cut placed on this quantity in the 
analysis [82] 



X 

Figure 8.2 Flowchart of the iterative procedure for tuning the simulation to real 

data [831 

Figure 8.3 Number of tripped or disabled channels in the //Tr for the angular decay 

coefficient analysis [85l 

Figure 8.4 Number of tripped or disabled channels in the /ilD for angular decay 

coefficient analysis [851 

Figure 8.5 Comparison of total cluster charge in the North /xTr from simulation 

and data IHTI 

Figure 8.6 Comparison of total cluster charge in the South /xTr from simulation 

and data IHTI 

Figure 8.7 Comparison of peak strip charge over total cluster charge in the North 

/xTr from simulation and data [88l 

Figure 8.8 Comparison of peak strip charge over total cluster charge in the South 

^Tr from simulation and data [Ml 

Figure 8.9 Comparison of tracking residuals in the North /iTr from simulation and 

data [901 

Figure 8.10 Comparison of tracking residuals in the South ^uTr from simulation and 

data [901 

Figure 8.11 Average efficiency for pairs of larocci tubes in the /xID plotted against 

BBC rate [U 

Figure 8.12 Comparison of azimuthal track projections between simulation and data. WI\ 
Figure 8.13 Comparison of radial track projections between simulation and data. . [93l 
Figure 8.14 Total counts and combinatorial background for an example pT range 

and fits to the signal [Ml 

Figure 8.15 Acceptance x efficiency plotted against transverse momentum and ra- 
pidity [961 

Figure 8.16 Acceptance corrected J/ip yield plotted against transverse momentum 

and rapidity ^7\ 



XI 



Figure 8.17 Angular decay coefficients for J/V' — )• fJ''^fJ'~ [991 

Figure 8.18 Magnitude of the angular decay coefficients in four reference frames 

integrated across pT and xp IIOOI 

Figure 8.19 Frame- invariant angular decay coefficient, A, from measurements in four 

different reference frames IIUII 

Figure B.l Relative acceptance for cos'd-<p for various ranges in pT I12UI 

Figure B.2 Relative acceptance for cos'd-<p for various ranges in xp 11211 



CHAPTER 1. Overview 



The J /ip meson, a bound state of a charm and an anti-charm quark, was discovered in 1974, 
yet the mechanism for the binding of the two quarks is stih largely unknown. Several models 
have been proposed for the so-called production mechanism, but none of them satisfactorily 
describe all available data. Because the J /ip is composed of two heavy quarks, the production 
of the quarks should have a different energy scale from their relative motion, and production 
may be understood by separating the two scales using non-relativistic QCD (quantum chromo- 
dynamics). Such an understanding would shed light on one of the basic components of QCD: 
hadronization. In this document I will present two novel measurements sensitive to the J/ip 
production mechanism using the PHENIX detector. 

The first measurement utilizes polarized proton beams at RHIC for a transverse single spin 
asymmetry of J/V' mesons, sensitive to both the internal structure of the colliding protons 
and to the production mechanism of the J /ip. This is the first time such a measurement has 
been made, and the resulting asymmetry is inconsistent with zero, a fact which has strong 
implications for understanding of both the proton and the J /ip. 

The second measurement is a determination of the angular distribution of decay leptons 
in the J /ip rest frame. Measurements of such distributions have been made in the past, but 
very few have determined all relevant coefficients or considered the coordinate system used in 
the J /ip rest frame. Using various coordinate systems and measuring all relevant coefficients 
will provide a much richer understanding of the J /ip at a phenomenological level and should 
provide guidance to theorists in pursuit of the production mechanism. 

The organization of this document will be as follows: In Chapter [2| I will give a short 
review of QCD, both to give the reader a bit of background on the subject and to introduce 



terminology which will be used throughout the document. Chapter |3] will discuss the structure 
of unpolarized protons, and Chapter |4] will present the current theoretical and experimental 
understanding of transversely polarized protons. In Chapter [5} I will give an overview of the 
J/ip production mechanism as it is currently understood. Chapterp^will present a description of 
the RHIC complex and the detectors making up the PHENIX experiment, detectors which are 
used in the analyses discussed in subsequent chapters. In Chapter [7) I will discuss the methods 
for measuring the transverse single spin asymmetry of J/ij: mesons, present the results of such 
a measurement, and give a short discussion of the results. Chapter [8] will present a similar 
discussion of the angular decay coefficients for J /ij) mesons. Finally, Chapter ^ will present 
conclusions which can be drawn from the two measurements and a possible path to an even 
richer understanding. 



CHAPTER 2. Quantum Chromodynamics 



The prediction and discovery of the J /ip meson and other bound states of charm and 
anti-charm quarks are the culmination of a vast theoretical foundation: the Standard Model 
of Particle Physics, Quantum Chromodynamics (QCD), and the parton model. In order to 
understand the relevance of the J/tp meson in enriching our understanding of QCD, we need 
to get our footing in the theoretical framework which predicted its discovery. 

2.1 The Standard Model 

The Standard Model of Particle Physics is arguably the most successful theoretical model 
in contemporary physics. It was first proposed by Glashow (95) as an attempt to describe the 
hundreds of new particles being observed at moderate and high-energy particle colliders in the 
1950s and 1960s and was later modified and improved by Weinberg and Salam (150; 144). The 
model is based on the idea that the majority of the matter we observe is composed of a relatively 
small number of common constituents called leptons, quarks, and gauge Bosons (93; 88). 
The Standard Model describes the interactions between the constituent particles in a single 
elegant expression, the Standard Model Lagrangian, derived from fundamental symmetries 
using group theory, and it has provided some of the most precise predictions in the history 
of theoretical physics. While it is not in the scope of this document to present a detailed 
discussion of the predictions and experimental status of the Standard Model (I will instead 
provide a brief introduction), a comprehensive guide can be found in the standard reference of 
particle physics (31). 

The integer spin mediators of fundamental interactions in the Standard Model are known 



as gauge Bosons (Table 2.1), and the fundamental spin-^ particles are the leptons and quarks 



(Table 2.2). Particles composed of two and three quarks are named mesons and baryons 



respectively, collectively known as hadrons. 



Table 2.1 Properties of the Standard Model gauge Bosons. Dependence 



of the strong force on distance will be discussed in Section 2.3 
From (31). 



Interaction Mediated 


Electro- Weak 


Strong 


Electromagnetic 


Weak 


gauge Boson 


7 (Photon) 


w+i- 


ZO 


gluon 


mass (GeV/c2) 





80.398±0.025 


91.1876±0.0021 





Electric Charge 





+1/-1 








Strength Relative to EM 10-^*m 
coupling for two up quarks at 3 X 10— ni 


1 


0.8 


0.8 


25 


1 


10-4 


lo--* 


60 



Table 2.2 Properties of the Standard Model quarks and leptons. 
From (31). 



Leptons 


Quarks 


Flavor 


mass (GeV/c2) 


Electric Charge 


Flavor 


mass (GeV/c2) 


Electric Charge 


electron 
neutrino 


< 7 X 10-9 





up 


.... ^ , .0-3 


2 
3 


electron 


0.54858 X 10-^ 


-1 


down 


.04 !J;- X .0-3 


1 
3 


muon 
neutrino 


< 0.0003 





charm 


1 27 +"■'" 

' -o.u 


2 
3 


muon 


0.10566 X 10-3 


-1 


strange 


105 % X 10-3 


1 
3 


tau 
neutrino 


< 0.03 





top 


m. !;;J 


2 
3 


tau 


1.77684 ±0.17 


-1 


bottom 


-° ^ 


1 
3 



The overriding symmetries of the Standard Model are given by the group SU(3)(S'SU(2)(8)U(1) 
The SU(3) group provides the symmetries for the strong interaction of quantum chromody- 
namics (QCD), SU(2) for the weak interaction, and U(l) for the electromagnetic interaction of 
quantum electrodynamics (QED). The number of gauge Bosons for a given interaction corre- 
sponds to the number of dimensions of the gauge group: A^^ — 1 for SU(A^). This means that 
there is one gauge Boson (the photon) for the electromagnetic interaction, 3 gauge Bosons (the 
W^ ^ W~ , and Z^) for the weak interaction, and 8 gauge Bosons (8 types of gluons) for the 



strong interaction. 

2.2 Quarks and Gluons 

In QED, stable particles have only three possible charges: positive, negative, or zero. In 
QCD, however, there are 6 charges and 6 anti-charges, and all stable particles which can be 
observed in isolation have a net charge of zero. The charges are called red, blue, and green with 
anti-charges anti-red, anti-blue, and anti-greenjj They are collectively referred to as 'colors' 
and 'anti-colors,' and stable particles are said to be 'colorless,' while charged particles are 
'colored.' All gluons in QCD carry an equal number of colors and anti-colors, and it is these 
charges which determine the type of gluon. The eight gluon types are: 

rf rb rg ' 
g« = Aj If bb bg 

\ gr gb gg J 

for i = 1,2, ...8, where the Aj are the Gell-Mann matrices, which are generators of the SU(3) 
group, and r, b, and g are red, blue, and green respectively (the bar denotes anti-color). The 
eight color configurations for gluons are collectively known as color-octet states. All gluons are 
colored so that none inhabit the color-singlet statqj 

— {rr + bb + gg). 

Quarks carry a single color (r, b, or g), while anti-quarks carry a single anti-color (f, b, or g). 
Leptons are colorless and do not experience the strong interaction. 

The fact that gluons carry color charge makes the strong interaction considerably more 
complicated than electromagnetism, where the photon has zero charge. An interaction wherein 
the gauge Bosons have zero charge, like electromagnetism, is called Abelian. An interaction 
wherein the gauge Bosons are charged, like the strong or weak interactions, is called non- 

Abelian. The non-Abelian nature of QCD implies that gluons interact not only with quarks 

^To have a net color charge of zero, a particle can either have an equal number of colors and anti-colors, like 
red and anti-red, or it can have one of each color or anti-color, red, blue, and green or anti-red, anti-blue, and 
anti-green. 

^The definitions of a color-octet and color-singlet state will be particularly useful when we talk about J/ip 
meson production in Chapter pi 



but also with other gluons. The interaction between gluons leads to a theory wherein quarks 
and gluons cannot be found in isolation but are confined to bound states (confinement) by a 
force which grows stronger as they move away from other quarks and gluons and weaker as 
they approach one another (asymptotic freedom). 

2.3 Asymptotic Freedom and Confinement 

The coupling constants of the Standard Model determine the strength of a given interaction. 
They are not 'constants,' however, as they depend on the mass of the exchanged gauge Boson. 
If a gauge Boson is short-lived (virtual), the Heisenberg uncertainty principle allows for its 
mass to fluctuate away from the nominal value, and even a virtual photon or gluon is allowed 
to have some mass. The mass of the gauge Boson, real or virtual, sets the energy scale Q of 
the interaction, defined as Q^ = —(mass) . 

If we know the QED coupling aem at a scale Q^ = m?^ the coupling can be calculated at 
a much larger scale q^ ^ m? by 

aem{q^) = a,^{m')[l + "^^""'^ log 4 + 0(«L)]- (2.1) 

For QCD we have 

as{q^) = as{rr?)[l + ^^^ log 4(2^/ " H^) + ^(4)], (2-2) 

A-Zil\ III 

where n/ is the number of (active) quark flavors, which depends on the energy scaler! and 
N is the number of colors. There are 3 colors and no more than 6 active flavors so that as 
decreases with increasing Q^, while aem increases. To put things a bit more simply: Q is 



he 



inversely proportional to distance d ~ tj, and Eq. 2.2 implies that the strong interaction is 
weaker for two quarks which are close together than it is for two quarks which are far apart. 
This surprising discovery, known as asymptotic freedom (139; 101), leads to the conflnement 
of quarks into bound state hadrons. 

A non-relativistic approximation for the interaction potential between a quark and anti- 
quark can be written (104) in terms of r, the radial distance between quark and anti-quark. 



^ Since the quarks have very different masses, a flavor of quarks is only relevant at an energy scale greater 
than or approximately equal to the rest mass energy for that flavor of quark. 



as 



V,,{r) = -^ (l - ^ + g(33 - 2n,)[ln(H + 1e]) + ^ (2.3) 

where as is the strong couphng constant, Uf is the number of active flavors, eo is the permit- 
tivity of free space, ^ and 'Je depend on the renormahzation schemer The color factor / is , 
for a color-singlet configuration and — i for a color-octet configuration. The divergence of the 
potential in Eg. |2.3| for r — )• oo, as well as the (related) increase of as at large distances, mean 
that colored objects (quarks, gluons, and colored hadrons) cannot be found in isolation. The 
difference in sign between the color factors for the color-singlet and color-octet configurations 
implies that the quarks making up a stable hadron can only be in a color-singlet configuration. 
The fact that quarks and gluons must be confined to colorless hadrons in QCD implies that a 
single quark or gluon can never leave a collision in isolation. Instead, it is energetically favorable 
for a parton leaving a collision to continue pulling colorless pairs of quarks and gluons from the 
vacuum until a totally colorless state can be reached. This process is known as fragmentation, 
and the binding of partons into colorless objects is called hadronization. Fragmentation cannot 
be calculated from first principles in QCD. Instead, the properties of fragmenting quarks 
and gluons are parameterized by fragmentation functions determined from experimental data. 
The parameterized measurements use groups of particles, called jets, observed in a particle 
detector every time a quark has fragmented, and the parameterization depends heavily on the 
experimental definition of a jet. 



*The theory of QCD contains several divergences, and a renormalization scheme is a method to handle those 
divergences in a rigorous way. Renomalization in QCD was first introduced in the 1970s by 't Hooft (149) 



CHAPTER 3. The Proton 



The particle collider is one of the most common tools for studying elementary particles. 
Collider experiments measure the flavor of outgoing particles and the directions in which they 
are produced from a collision. These measurements can then be compared to theoretical 
predictions in order to develop a combined understanding. 

Many contemporary experiments study particle production by colliding two beams of pro- 
tons. As we will show, the proton is a composite object, made of quarks and gluons. To develop 
an understanding of the particles produced in proton collisions, like the J/V' meson, we need to 
have a thorough understanding of the structure of the proton itself. Otherwise, quantities like 
cross-sections from proton collisions cannot be predicted, leaving little theoretical guidance for 
experiment. 

3.1 Cross-Sections and Luminosities 

In scattering experiments, it is necessary to quantify the number of particles measured 
so that results can be compared between separate experiments. This is done in part with 
differential cross-sections ^ where dQ is a unit of solid angle and 

is the total cross section. The cross-section depends on the luminosity £ as cr = ^^, and for 
a beam of particles in a circular collider 

C = fn^j- (3.2) 



where / is the crossing frequency and n the number of bunches of colhding particles in one 
beamjj A'^i is the number of particles per bunch in one beam, while A'^2 is the number of 
particles per bunch in the other beam, and A is the cross-sectional colliding area of the beams. 
The differential cross-section can be measured with respect to any physical quantity, but 
one of the most basic is the invariant cross section with respect to the energy E and three- 
momentum p. Through a change of variables, this quantity can be rewritten as 

j-.d'^cr 1 (Pa 

dp^ 2-7TPT dydpT 

where y is the rapidity and px the transverse momentum of the particle. 

For measurements of particle production, rapidity is preferable to 0, the angle relative to 
the beam axis, because it is additive while maintaining Lorentz invariance. Rapidity is defined 
as 

where E is the total energy of the particle and pz its momentum along the beam axis. It is 
often convenient to use the pseudo-rapidity, a quantity which is equivalent to rapidity if the 
particle is massless but depends only on the outgoing angle of the particle: 

1 1 (\P\+Pz 

rj = - m 



2 \\p\-Pz 



= -In (tan- J. (3.5) 

Experimentally, both the luminosity and the detector efficiency e must be known in order to 
determine a cross section. For example, the differential cross-section as a function of transverse 
momentum pT can be written as 

#^ = 7^^ (3.6) 

apT Ce dpT 

where dN is the number of particles measured in a transverse momentum range dpT, and C is 
the provided luminosity. 



^Beams of particles in circular accelerators are not continuous but are composed of a number of equally 
spaced 'bunches.' 



10 
3.2 Deep-Inelastic Scattering 

Scattering of leptons on protons has long provided precise information on the substructure 
of the proton. The charged point-hke structure of leptons makes them powerful tools for 
probing the electronic structure of the much larger proton. Early measurements compared the 
total scattering cross-section for ep — )• eX with calculations from several ansatz for a point-like 
charged proton structure: 

• The 'Mott cross-section' for scattering of a spin-1/2 electron on a classical potential 

• The 'Dirac cross-section' for scattering of a spin-1/2 electron on a spin-1/2 point charge 

• A modified Dirac cross-section, which takes into account the anomalous proton magnetic 
moment. 



The measured cross-section does not agree with any of these point-like ansatz (see Fig. 3.1), 
implying that the proton is not a point-like objectjj 

The quantum numbers of the proton (charge and isospin) indicate that its dominant con- 
stituents in the quark model ought to be two up quarks and one down quark, called the valence 
quarks. While this is a good approximation for collisions which are barely inelastic, high en- 
ergy collisions begin to probe the QCD vacuum structure of the proton, and a 'quark sea,' 
consisting of all possible flavors of quark and anti-quark, is observed along with a very large 
density of gluons. Since the details of this structure cannot be calculated from first principles, 
the probability distributions for finding a gluon or a specific flavor of quark or anti-quark are 
determined experimentally with deep-inelastic scattering (DIS). 

In DIS, a lepton is scattered on a proton with enough energy to probe its internal structure, 
and the outgoing lepton is observedr The scattering resolution is ~7y where Q^ = —q^q^ and 
q = k' — k with k the 4-momenta of the incoming lepton and k' the 4-momenta of the same 



^Admittedly, the anomalous magnetic moment of the proton already provides evidence that it is not a point- 
like particle, but the verification of this conclusion with electron scattering was quite important historically. 

^Semi-inclusive deep-inelastic scattering (SIDIS) is equivalent to DIS except that other outgoing particles 
are measured in addition to the scattered lepton. 
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Figure 3.1 Scattering of a 188 MeV electron beam on a hydrogen target 
showing that the proton is a composite object. From the body of 
work that won Hofstadter and Schiff the 1961 Nobel Prize (108). 



lepton after scattering (Fig. 3.2). Because the proton is approximately 1 Fermi in diameter, 
quarks become the relevant degrees of freedom at (5~0.5 GeV. 

The scattering cross-section of a lepton on a proton with 4- momentum P is da oc L^yW^'^ 
where 

L^,v oc y^ {u{k\ \')j^u{k, A)) {u{k', X')^t,{0)u{k, A))* , 



W^^" 



oc 



A, A' 



E 



^.i.. 



27r 



p,c7\[Jiio,Jum\p, 



a 



are the leptonic and hadronic tensors respectively, and J is a matrix describing the effective 

vertex between the virtual photon and the protonrl The helicities (scalar product between the 
''An effective vertex uses a single coupling to describe a much more complicated interaction. 
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P,M- 

Figure 3.2 Illustration of commonly used variables in DIS. 

spin and momentum vector) of the incoming and outgoing lepton are A and A respectively, 
while a is the helicity of the incoming proton. W^i, is considered universal, meaning that the 
same W^u can be applied for scattering of a proton on any particle, not just a lepton. The 
total scattering cross-section for a proton-proton collision, for example, is da oc Wf^yW^^ . 

3.3 Structure Functions and Bjorken Scaling 

For unpolarized DIS using electrons or muons, the hadronic tensor can be parameterized 
as 

W^, = g^,,Wi{xB, Q") + P^PuW2{xB, 0") (3.7) 

where P^ = {P^ '^q^)/M, xb = ^m o ' ^ ^^ ^'^^ mass of the proton, and Ij the inverse metric 

tensor. 

The quantity xb (called the Bjorken-x) is especially convenient, because for large proton 
momenta it can be interpreted as the fraction of the hadron momentum carried by a quark. 
An additional convenience is that at large energy transfer g*^ — t- cxd and fixed x^, the structure 
functions W\[xb^ Q'^) and W2{xb, Q'^) depend only on xb and not Q^ so that MWi{xb, Q^) = 
Fi{xb) and ^W2ixB,Q^) = MF2{xb) (55). 

Fi{xb) and F2{xb) are called the 'unpolarized proton structure functions,' and they can 
be related to the probability for the scattering of a polarized virtual photon on a quark in 
the hadron. In the parton model, the longitudinal structure function Fl = F2 — 2xbFi (cor- 
responding to a longitudinally polarized photon) can be written as Fi = -r-^ — cfl{xb^Q'^)j 

47r ttem 

where gl is the cross section for interaction between a quark and a longitudinally polarized 
photon. Similarly, the transverse structure function Ft = xbFi (corresponding to a trans- 
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versely polarized photon) can be written as Ft = -r-^ — o'tIxbjQ'^), where ar is the cross 

section for interaction between a quark and a transversely polarized photon. Callan and Gross 

have shown that <Ti/o"T scale as 1/Q^, meaning that F/, — )• as g" — t- oo (65) so that Fi and 

F2 are related at large momentum transfer. 

In proton-proton collisions, the quantities xi and X2 are defined (138) to be approximately 

equivalent to xb at leading order in as as 

Af „ M „ , , 

xi = -^ey, X2 = -^e-y 3.8) 

for partons from the proton travelling in the positive and negative z direction respectively, 
where M is the particle's mass, y its rapidity, and ^/s the center of mass collision energy. If 
the measured particle is sufficiently massless, M can be replaced with the transverse momentum 
Pt and y with the pseudo-rapidity r]. 

The parameterization of Eq. |3.7| along with Bjorken's observation that the structure func- 
tions do not depend on Q^ as g*^ — t- 00 are very powerful due to the fact that they can be easily 
interpreted in term of quarks. At significantly large energies in electron-proton collisionsjjthe 
electron probes the structure function F2 of a proton as 

F2ix)=x"^e^[q{x)+q{x)], (3.9) 

where q{x) {q{x)) is the probability to find a quark (anti-quark) of flavor q at momentum 
fraction x. Similar expressions can be determined for DIS using other lepton probes (i.e. 
uP — ^ vX and e'^P — t- e^X), and these measurements together determine both q{x) and q{x). 

3.4 Altarelli-Parisi Evolution 



Part of the reason that Eq. 3.9 and the Callan-Gross relation work as q^^-oo is that 
gluons have very little effect on high energy quarksjj At lower q^, gluons become much more 

important, and Bjorken scaling no longer applies. The contribution of gluons to the structure 

^The actual approximation is that the perturbative expansion is taken only to leading order in Qem, the 
electromagnetic coupling constant. 

^Another reason is that the momenta of most quarks is preferentially along the direction of the proton's 
motion, meaning that the effects of transverse motion in the proton are negligible. Transverse partonic motion 
inside the proton will be discussed at length in Chapter Wl 
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function F2 follows Altarelli-Parisi evolution, a set of differential equations relating the quark 
distributions to gluon distributions: 

dg{x,Q) as(Q^) 



d log Q vr 



/' -{Pg^z) J^fe(^, Q) + qi^ Q)] + P,^g{z)gi^Q)}, 

Jx Z Z Z Z 



where the g{x,Q), q{x,Q), and q(x,Q) denote the probability for the scattering of a gluon, 
quark, and anti-quark respectively (30). The splitting functions Pj-t-i are 



^5^.(^) 



4 / l + (l-z)2 
3 V z 

^ '2 I ^1 ^^2^ 



P^^gi^) = ^{z' + {l-z 

where n/ is the number of (active) quark flavors, and [F{z)]^ is defined as 

dzfiz)[F{z)]+= I dz{f{z)-f{l))F{z). 
Jo 

The Altarelli-Parisi equations can be used along with the measured q{x) and q{x) distributions 



from Eq. 3.9 to determine both the quark and gluon parton distribution functions (PDFs). 

The structure function F2 has been extensively measured, and a compilation of results 
can be found in Fig. |3.3[ In order to determine the PDFs from these measured structure 
functions, a global analysis must be performed, taking into account information from several 
scattering processes at various energies. Global analyses are continually updated in order to 
include the most recent structure function data and to improve on methods for the parton 
model fits (140; 51; 130). An example of a recent fit from the MSTW collaboration (130) is 
shown in Fig. |3.4| where it is clear that the composition of the proton is well determined for a 
large range in x. 
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Figure 3.3 Compilation of measurements of the F2 structure function 
(from (31)). Deviations from the flat line (which is indicative of 
Bjorken-scaling) are due to Altarelli-Parisi evolution at small-x 
and elastic scattering contributions at large-x. 
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MSTW 2008 NNLO PDFs (68% C.L.) 
Q2 = 10GeV2 




Figure 3.4 Parton Distribution Functions from the MSTW collabora- 
tion calculated at next-to-next-to-leading order in Os for 
q2 = 10 GeV2 on the left and Q^ = 10^ GeV^ on the right 
(from (130)). 
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CHAPTER 4. Transversely Polarized Protons 



The discussions of Chapter |3] focused on unpolarized protons, but an even richer under- 
standing of the proton structure can be developed by studying colhsions of polarized protons. 
The structure of a transversely polarized proton is different from that of a longitudinally 
polarized proton, because partons are highly relativistic so that their boosts and rotations 
do not commute. In this work, I will only discuss the structure of a proton with spin ori- 
ented transverse to its momentum (transverse polarization), but a great deal of understanding 
has also been developed from the study of protons polarized along their momentum direc- 
tion (96; 97; 78; 79). 

Early experiments with transversely polarized protons yielded surprising results, and the 
theoretical understanding continues to develop. Before we discuss the measurements and the 
associated theory in detail, it is necessary to define the basic observable: a transverse single 
spin asymmetry (SSA). Put simply, transverse SSAs are an imbalance of the production cross- 
section to one side of the proton spin relative to the other. 

For proton-proton collisions in a collider, the asymmetry is measured using particle yields 
in the left and right hemispheres from opposing spin orientations of a single polarized beam. 
The term 'left' is defined as the axis which forms a right-handed coordinate system between 
the polarized beam direction and one of the proton spin orientations, denoted as t- If the 
momentum vector of the outgoing particle is p and the colliding proton momentum is P with 
spin 5, left is defined as p • ( 5 x P j > 0, and the transverse SSA on the left is 

^^ = |4^ (4.1) 

where a\ ((T|^) represents the production cross-section in the left hemisphere with beam polar- 
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ized in the t (i) direction, and V is the beam polarization. An overall minus sign is required 
for A]^ on the right. 

The geometric scale factor / takes into account the convolution of an azimuthal asymmetry 
with an incomplete detector acceptance. We assume that the modulation of A^ is sinusoidal 
in azimuth (a shape which will be justified in the following sections) and ensure that we are 
measuring the amplitude of such a modulation by defining 

i r ■ \~^ 

/ f((ii) .sin cbdfh 



/ = 



e{(j)) sin (j)d(j) 



(4.2) 



for An on the left, where €{(j)) is the efficiency for measuring a particle with azimuthal angle 
(j) between p and S. The limits of integration correspond to the hemisphere in which the 
measurement is made. 



Measurements of A^ (Eq. 4.1 ) for various event topographies currently provide most of the 



basis for understanding the structure of a transversely polarized proton. 

4.1 Experimental Observations 



QCD naively predicts that An should scale as rnqj ^fs^ where niq is the mass of the scattered 
quark from the polarized proton (117). Because a large fraction of the proton is composed 
of light up and down quarks, the asymmetry should be 0(10"'') for -y/s =20 GeV collisions, 
but experimentally measured asymmetries have consistently been 0(10"^). This contradiction 
has led to a great deal of theoretical and further experimental activity, and the theoretical 
understanding of large measured asymmetries is still far from complete. 

Unexpectedly large transverse SSAs of up to ~ 40% were first observed in pion production 

from a ]r bearq^ incident on a liquid hydrogen target in 1976 at the Argonne zero- gradient 

synchrotron (ZGS) (120) and subsequently in a number of experiments using hadronic collisions 

over a range of energies extending up to y^ = 200 GeV. Several examples of such measurements 

are: 

'^The t denotes a transversely polarized proton 
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tt"*"' production from a 22 GeV/c polarized proton beam incident on C at the Brookhaven 
AGS (29) 

vr*^ production from a 24 GeV/c proton beam incident on a p^ target at the CERN 
experiment PS 141 (38) 

71-+/- production from a 200 GeV/c polarized proton beam incident on a liquid Hydrogen 



target at the E704 experiment (5) (Fig. 4.1) 



vr'^ production from a 200 GeV/c p beam incident on a p' target at the Fermilab E581 
experiment (6) 

7r+' and K^' production from p^+p collisions with ^/s =62.4 GeV at the BRAHMS 



experiment at BNL (42) (Fig. 4.2) 



vr^ production from p'+p collisions with ^/s =200 GeV at the STAR experiment at 



BNL (7; 2) (Fig. 4.3) 



vr^ production from p'+p collisions with y^ =200 GeV and ^/s =62.4 GeV at the 



PHENIX experiment at BNL (18; 22) (Fig. 4.5) 



The asymmetries have also been observed in SIDIS using p^ (23; 81) and polarized deuteron 
(27; 127; 26) targets. 

Transverse SSAs are typically measured as a function of Feynman-x, xp = xi — X2, where 
xi is the X of the scattering parton from the polarized and X2 from the unpolarized proton. 
The asymmetries become non-zero around xf=0.2 and are zero for xp < 0, statements that 



hold true for collisions at ^^=19.4 GeV, (Fig. 4.1), ^^=62.4 GeV (Fig. 4.2), and ^=200 GeV 



(Fig. 4.3). Because the up and down valence quarks dominate for x > 0.1, we can infer that 
the transverse SSA is dominated by effects from valence quarks in the polarized proton. 

The theoretical explanations of large transverse SSAs typically suggest additional parton 
distribution functions which describe correlations between the polarization of the proton and 
the partons which make up that proton: 

The Sivers Function. Distribution of unpolarized quarks in a transversely polarized proton. 
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Figure 4.1 Fermilab E704 transverse SSAs for n^, n^ and vr from p+p col- 
lisions at y/s = 20 GeV as a function of Feynman-x (from (5)). 
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Figure 4.2 Brahms transverse SSAs for vr"'' and vr^ from p+p collisions at 
-y/s = 62.4 GeV as a function of Feynman-x for various ranges 
of transverse momenta (from (42)). 



The Boer-Mulders Function. Distribution of transversely polarized quarks in an unpolar- 
ized proton. 

Transversity. Distribution of transversely polarized quarks in a transversely polarized proton. 

In addition to these distributions, transverse SSAs can be produced by the Collins effect (the 
correlation between the spin of a polarized quark and hadrons in a fragmenting jet). 

While these distributions each describe some aspect of the hadronic structure, the mecha- 
nisms for creating asymmetries are quite different. In reality, it is likely that no single effect 
causes the transverse SSA but that effects from all of these additional parton distribution 
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Figure 4.3 STAR transverse SSAs for vr'' from p+p collisions at 
^/s = 200 GeV as a function of Feynman-x (from (2)). 



functions contribute. 



4.2 The Sivers Effect 



In 1989, Dennis Sivers proposed that partonic motion transverse to the proton's momen- 
tum, kT, ought to be taken into account when calculating PDFs for transversely polarized 
protons (145; 146). He suggested that this could be done by creating transverse momentum 
dependent parton distribution functions fgipt{x, kx, St) (TMDs) as opposed to the usual f1{x) 
used for unpolarized protons, where q corresponds to the flavor of quark, anti-quark, or gluon. 
The new TMD includes a so-called Sivers function, A^/g/pt(x, /ct): 

fg/pt{x, kr, 5t) = fUx, kr) + ^A^/g/pt(x, kT)ST ■ (P x fcr) (4.3) 

for a parton with transverse momentum kx and spin component St transverse to the three- 
momentum P of the proton. Alternatively, the TMD can be written as 

1 



fg/ptix, kr, St) = f{{x, kr) + -A fq/ptix, kr) sin((/>5 - ip) 



(4.4) 
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where (ps is the azimuthal angle of the proton spin and ip the angle of the partonic kx- Because 
the Sivers function is weighted by a sinusoidal modulation, it can be observed as a difference be- 



tween TMDs, corresponding to the numerator of Eq. 4.1 , or directly as a sinusoidal modulation 
in particle production. 

The Sivers effect is not universal. Instead, it should have the same magnitude and opposite 
sign in SIDIS relative to Drell-Yan production (quark/anti-quark annihilation into a pair of 
leptons) from p+p collisions (75). There has not yet been an experimental verification of this 
sign-change, but such a verification will be quite important in proving the existence of a Sivers 
effect. 



Unfortunately, Eq. 4.3 and the associated formalism is not valid for hadron production 
from hadron-hadron collisions. It was shown in 2007 that factorization is violated for the 
TMD approach to the Sivers function in the case of back-to-back hadron production from 
hadron-hadron collisions (73). Furthermore, in 2010, such violations were shown to occur for 
any hadronic production from hadron-hadron collisions (143). 

An alternative approach to calculating the transverse SSA caused by unpolarized quarks 
in a transversely polarized hadron was developed by Qiu and Sterman, wherein the transverse 
SSA is generated by an incoming or outgoing parton which exchanges a gluon with a parton 
in the initial state hadron remnant (141). Kanazawa and Koike found that the same effect 
could occur for the exchange of a gluon with a parton in the final state remnant (116). The 
collinear approach, which does not violate factorization, orders the perturbation theory in 
powers of ( -jy^ ] \h] i where P is the four momentum of the polarized hadron, Q^ the 



momentum transfer, q the momentum of the quark in the polarized hadron, d the dimension 
and s the spin of the operator in a given term of the expansion. The size of the contributions 
decrease with the twist, d — s, of the contributing operators, and it is the twist-3 terms, called 
multi-parton correlation functions, which cause the transverse SSA. For processes involving 
only gluons we will refer to these twist-3 terms as trigluon correlation functions. 

Because the collinear approach expands about i, the perturbation theory is valid for 



23 

collisions with transverse momentum transfer q± and Q ^ Aqc"/) PI while the TMD approach 
requires a separation between scales of the transverse and longitudinal momentum transfer so 
that q± <^ Q. The two formalisms have been shown to be equivalent in their region of overlap 
Aqcd < g± « Q (113)0 

4.3 Trans versity 

The transversity distribution describes the kinematics of transversely polarized quarks in 
a transversely polarized hadron (142; 112). Transversity cannot be responsible on its own for 



the large transverse SSAs discussed in Section 4.1 , because quarks from an unpolarized hadron 



in hadron-hadron collisions will naively scatter with an equal probability in all polarizations. 
Likewise, in SIDIS, the virtual photon interacting with quarks in the polarized hadron will not 
have a preferred polarization direction. Instead, the transversity distribution must couple with 
another distribution to create a transverse SSA. The distribution coupled with transversity can 
either be a property of the initial state of the unpolarized hadron, such as the Boer-Mulders 



function (Section 4.4), or it can be a final state distribution, such as the Collins function 



(Section 4.5). 



Soffer has shown that the transversity distribution is bounded by measurements of the 
unpolarized and longitudinally polarized PDFs (147). The model-independent 'Soffer bound' 
requires that 

|Atq(x)| < -{q{x) + Aq{x)) 

where Axqix) is the transversity distribution for a quark of flavor q, q{x) is the unpolarized 
PDF, and Aq(x) the longitudinally polarized PDF. Because gluons are massless, they must 
have a definite helicity, which means that the transversity distribution for gluons must be zero. 



^J^QCD depends on the renormalization scheme and does not have a clear physical interpretation, but it is 
roughly the energy scale at which quarks bind into hadrons. 

^Throughout the rest of the document we will use the term 'Sivers effect' when discussing the TMD and 
colhnear approaches collectively. This should not be confused with term 'Sivers function,' which is only appli- 
cable to the TMD approach. 
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4.4 The Boer-Mulders Effect 

Expanding on the idea of the Sivers effect, Boer and Mulders suggested that transverse 
SSAs could also be generated by polarized quarks in an unpolarized hadron (57). The Boer- 
Mulders function, h^ Jx^kx), corresponds to a modification of the unpolarized PDFs which 
takes into account the polarization of the constituent quarks in an unpolarized hadron. The 
parton distribution function for a polarized quark in an unpolarized hadron becomes 



1 / ST-[P><kT 
/,t/p(x, Tt, St) = -\ ffix, kr) - hi^ix, kr) \^ '- | (4.5) 

where M is the mass of the hadron (49). A non-zero Boer-Mulders function in an unpolarized 
hadron coupled with transversity in a polarized hadron could potentially lead to transverse 
SSAs. Such a distribution should also create modulations in the azimuthal distribution of 
particle yields from unpolarized hadron collisions. 

4.5 The ColHns Effect 

In an alternative explanation of large transverse SSAs, Collins suggested that a transversely 
polarized quark should fragment differently than an unpolarized quark, causing an azimuthal 
asymmetry about the jet axis (74). If an outgoing quark has spin Sq, the unpolarized fragmen- 
tation function D\{z,pr) is modified by a Collins fragmentation function A'^ Df^/^{z,pT) for 
a hadron h with mass Mh and transverse momentum px, carrying a fraction z of the quark's 
momentum kg: 

D^/gt{z,PT,Sq) = D\{z,pt) + 2\pt\A^ Dj,ig,{z,pT)S, ■ (pt x fc,) . (4.6) 

If there is a non-zero transversity distribution in the polarized hadron, a quark going into 
the hard scattering will have some polarization, and it is likely that an outgoing quark will 
maintain that polarization. The polarization of the quark causes its fragmentation function to 
be modified in such a way that there is a sinusoidal modulation in hadron yields about the jet 
axis. 
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The Collins function can not be larger than the unpolarized fragmentation function it 
modifies, a constraint called the positivity bound. 

4.6 PDF Extractions and Global Fits 

The kinematics of polarized SIDIS are such that the Sivers and Collins functions create 
different angular distributions and can be separated experimentally. Instead of the one dimen- 



sional distribution given in Eq. 4.1, SIDIS experiments take the target polarization direction 



into account and measure the two dimensional distribution 

for a hadron of type h. V is the target polarization, f and J, the spin direction, (j)s the azimuthal 
angle between the interaction plane and the target spin, and (f) the azimuthal angle between the 
interaction plane and the outgoing hadron momentum. The Sivers moment {sin{(p — 4's))ut 
and Collins moment {sin{4> + (J)s))ut '^^'^ then be extracted by fitting the two dimensional 
distribution with 

J^h \ — ly\ 

-^ = {sin{(f) - (t>s))uT sM(P - (t^s) + {sin{(l) + 4>s))uT -TTi \ / \. sin{<^ + 4>s), (4.8) 



3.2) and 



where y = -^A (with P, q^ and k defined in Section 

,, ,_y\ {l-y){l + R{x,y)) 
Mx.y) = y + 2Mxa 

^^ Ey 

with R{x, y) the ratio of the longitudinal to transverse structure functions discussed in Sec- 
tion [331 

The Sivers moment is sensitive only to the Sivers function and can be extracted directly. 
The Collins moment, on the other hand, is sensitive to a convolution of the Collins effect and 
transversity, and another measurement is necessary to extract both distributions. 

A global fit to extract the Sivers function (36) was performed on SIDIS data with a hydrogen 
target from HERMES (23; 81) and a deuteron target from COMPASS (27; 127; 26). The first 
moment in k^ of the Sivers function 

A^4'l>t(^) = / ^'^T^ A^/,/pt(x, kT) (4.9) 
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from the fit is shown in Fig. 4.4 The data favors a non-zero first moment in SIDIS which is 
positive for up quarks, negative and nearly identical in magnitude for down quarksj^ 
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Figure 4.4 Sivers function from a global fit to SIDIS data. The solid line 
is the central value of the fit from (36) and the dashed line 
from (34) (Figure from (36)). 



Because an electromagnetic probe is used, the gluon Sivers and trigluon correlation func- 
tions are not well constrained by SIDIS measurements. Measurements of the transverse SSA 
for vr'^s at PHENIX (18) seem to prefer a small gluon Sivers function at mid-rapidity (33), 



and recent preliminary results (Fig. 4.5) confirm a small transverse SSA out to large pT- A 
numerical constraint from the data has not been made, and such constraints are complicated 
by the process dependence of vr*' production. At low px, vr*^ production from p+p collisions at 
-y/s=200 GeV is dominated by gluon-gluon collisions, but quark-gluon collisions become more 
important at moderate pT, and quark-quark collisions become dominant at sufficiently large 



PT (148) (see Fig. 4.6). 



In order to extract the Collins and Transversity distributions from SIDIS, those measure- 
ments are taken together with measurements of the Collins distribution from BELLE (1). The 
BELLE detector is located at the KEKB accelerator, which collides unpolarized e^ and e^ , and 



''Recall that the Sivers function should have an opposite sign for hadron-hadron collisions. 
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Figure 4.5 Preliminary PHENIX transverse SSAs for vr*^ from p+p colli- 
sions at -y/s = 200 GeV as a function of pT- 



the Collins function appears as a modulation cos(0i + (j)2) in the azimuthal angles (pi and (p2 of 
outgoing hadrons from back-to-back jets with respect to the reaction plane. Because there are 
no hadrons in the initial state, the measurement from BELLE is only sensitive to the Collins 
function and not to transversity. This means that the BELLE data can be used in a global 
fit with data from HERMES and COMPASS to extract both the Collins and Transversity 
distributions (37). 



Transversity distributions from the global fit (Fig. 4.7) are non-zero and smaller than the 
Soffer bound and again appear to have an opposite sign for down quarks relative to up quarks. 

Collins distributions are extracted separately for favored fragmentation, where the observed 
hadron contains a valence quark in the flavor of the struck quark, and unfavored fragmenta- 
tion, where it does not. The Collins fragmentation functions extracted from the global fit 
(Figure 4.8) show an increase of the Collins eff'ect with respect to z as we would expect Pj 

The Boer-Mulders function is a property of unpolarized protons and can, in principle, be 

extracted from unpolarized SIDIS or Drell-Yan. Recent measurements at COMPASS (115; 61) 

^If a single hadron carries most of the momentum of a jet, the azimuthal distribution created by the quark 
polarization will not be diluted by extensive fragmentation. 
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Figure 4.6 Incoming partons from processes contributing to vr^ produc- 
tion at mid-rapidity from p+p collisions at ^/s=200 GeV 
(from (148)). 



and HERMES (94) of the cos 2<j) modulation in azimuthal angle (p of hadrons with respect to 
the production plane from unpolarized SIDIS were used in a proof of principle global fit to 
extract Boer-Mulders functions (52). 

One complication to the extraction of the Boer-Mulders function is the Cahn effect (63; 64), 
a cos 2(f) modulation in hadron production due to non-zero kx in the proton. Fortunately, the 
Cahn effect has a different kinematic dependence than the Boer-Mulders effect in pT, kx, and 
Q'^, and the two can be separated in the fit. While the available data do not allow a full 
extraction of the Boer-Mulders effect, a fit can be performed assuming that the Boer-Mulders 
function is proportional to the Sivers function extracted in (36). Results of a fit which makes 
such an assumption to extract the first moment in /c-p of the Boer-Mulders function 



hfi'Hx)^ fd'kT:^^hi^{x,k'r) 



(4.10) 



are shown in Fig. 4.9 The fit prefers a Boer-Mulders function which has approximately the 
same magnitude for up and down quarks, twice as large as the Sivers function for the up 
quark jj 

^The plots shown in this document are a bit confusing on this last point, because the Sivers function shown 
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Figure 4.7 Transversity distribution from a global fit to SIDIS and e'^e~ 
data. The blue line represents the Soffer bound. Light uncer- 
tainty bands are those from (35), while dark bands are from (37) 
(Figure from (37)). 



The global fits performed for the Sivers, transversity, Collins, and Boer-Mulders distri- 
butions all prefer magnitudes which are inconsistent with zero, implying that all of these 
distributions play some part in the large Ajsr seen by experiments. These fits are still in their 
early stages, however, and a great deal more data is required to sort out which effects dominate. 
The gluon Sivers effect, for instance, is still poorly constrained, and further measurements from 
hadron-hadron collisions, like the measurement presented in Chapter [7j will be necessary to 
fully explain the large measured transverse SSAs. 



in Fig. |4.4| is scaled differently than the Boer-Mulders function in Fig. |4.9| A thorough discussion of the notation 
used for the various transverse PDFs and the relationship between different notations is given in (49). 
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Figure 4.8 Collins function from a global fit to SIDIS and e+e^ data. The 
dark blue line represents the positivity constraint. Light uncer- 
tainty bands are those from (35), while dark bands are from (37) 
(Figure from (37)). 





Figure 4.9 Boer-Mulders function from a proof-of-principle global fit to 
SIDIS data assuming that the Boer-Mulders function is propor- 
tional to the Sivers function from the global fit in (36) (Figure 
from (52)). 
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CHAPTER 5. The J/^ Meson 



5.1 Charmonium 

Analogous to the term 'positronium' (a bound state of e^e~), the term 'charmonium' refers 
to any bound state of a charm and anti-charm quark. The various charmonium states differ in 
their total angular momentum, charge conjugation, parity, and principal quantum number. It 
was proposed soon after the discovery of such bound states that a spectroscopic arrangement of 



the various resonances would be beneficial (80). Such an arrangement can be seen in Fig. 5.1 
and continues to be used in lattice QCD to successfully predict the masses of and transitions 
between states (77; 82). 

The most copiously produced charmonium state is the J/V'i a spin-1, parity odd, charge-0 
meson in an s-wave orbital state. The J /ij) meson was simultaneously discovered at BNL (46) 
and SLAC (47), where it was named the 'J' and '-0' respectively. A composite name has 
remained, although early papers differ on the order. Since the discovery of the J/^, many 
other charmonium resonances have been confirmed experimentally, a number of which can be 



seen in Fig. 5.2 It is important to note that many of the J /if} mesons observed at colliders 
are not directly produced from collisions but are the result of decays from other states in the 
charmonium spectroscopy or from decays of B mesons or heavier charmonium states jj Recent 
models estimate that approximately 30ibl0% of J /ij) mesons produced at RHIC come from Xc 
decays, and overall only 59ibl0% of the J /ifj mesons are directly produced (62). 



'^B mesons are composed of a bottom or anti-bottom quark along with an up, down, charm, or strange quark 
or anti-quark. 
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Figure 5.1 Various charmonium resonances and their decay channels 
(from (31)). Three Xc states, as well as the ip{2S), contribute 
to J/'ip production in addition to those directly produced. 



5.2 The OZI Rule 



Because gluons are colored and hadrons are not, interactions with a single gluon connecting 
the initial and final states ought to be suppressed. The OZI rule, posited by Zweig (153) and 
elaborated by lizuka (109) and Okubo (136; 137), generalizes this principle to show that there 
is a supression of any interaction that can be split into separate initial and final states by 
cutting through gluon lines. This simple observation has striking consequences for decaying 
bound states like charmonium, as it explains the suppresion of decays to hadrons relative to 
leptons and the associated narrow decay width (93.2ib2.1 keV (31)) of the J/V'- While we 
might expect diagrams for charmonium decays to hadrons like that in Fig. |5.3|to make up the 



entire decay width, the observed decays to leptons make up a significant fraction (Table 5.1). 
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Figure 5.2 Various charmonium resonances confirmed in data from several 
experiments. Tlie y-axis is the R value, defined as the ratio of 
the cross section of oppositely charged hadrons from a colliding 
e"*" and e~ at center of mass energy ^/s to the cross section for 
producing a pair of oppositely charged muons (figure from (31)). 




Figure 5.3 Example of a decay process suppressed by the OZI rule. 



5.3 Angular Distributions of Decay Products 

Charmonium systems are well described by hydrogen-like potentials having quantum num- 
bers n, I, and m corresponding to the principal quantum number, eignvalue for the total 
angular momentum J squared, and magnetic quantum number respectively. The J/ip is in an 
s-wave orbital state (/ = 0), meaning that it can't be 'polarized' in the sense that (J^) 7^ 0, 
because the quantum number m must be zerojj It is possible, though, for the spin to be 



■^The real-space wavefunction which solves a hydrogen-like Hamiltonian contains a radial piece mult iplied by 
Legendre polynomials, which are zero for \m\ > I. Likewise, the Wigner D-matrices -D^a in Eq. 



for \m\ > I. This is not coincidental, as the d-functions of Eq. |5.2| are related to the Legendre polynomials by 

4o = ^i(cos/3). 



5.1 
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Table 5.1 Dominant J/ip decay modes. Branching ratios are from (31). 



Decay Channel 


Percentage of all decays 


hadrons 


87.7±0.5% 


7* — ;> hadrons 


13.50±0.30% 


e+e" 


5.94±0.06% 


^+/i- 


5.93±0.06% 



preferentially aligned along some axis with an equal probability to be pointing either along 
or against that axis, and this alignment will cause a non-uniform angular distribituion for 
leptons from the decay. Here we present a derivation of the angular distribution of leptons 
from J /ip decays, borrowing heavily from the much more thorough treatment of angular decay 
distributions in (110). 

Suppose we have a particle 7 in eigenstates j and m of total angular momentum and 
magnetic sub-state respectively. The 7 decays into two particles a and (3 with momenta p and 
—p in the 7 rest frame and they have corresponding helicities A^ and A^ n The decayed state 
is characterized by a ket \p\a^p) where p is a unit vector in the direction of ^. This state can 
be expanded into angular momentum eigenstates of the parent 



IpA^A^) = Y^\jmKXp)^'^-^D^^^{^,§,i^) (5.1) 

where A = A^ — A^ and Z?^^(y3,'!9,'0) are the so-called 'Wigner /^-matrices' corresponding 
to matrix elements {jm'\R{Lp,'d,ip)\jm) of the rotation operator R{ip , •& , ifj) . The three Euler 
angles c^, ??, and V specify the direction of p, and the Wigner L>-matrices can be written in 
terms of the more convenient real- valued d-functions 

Dl^{^, ^, ^) = e~^"^^dl,{^)e-^'^, (5.2) 

some of which can be found in (31), as well as many other sources. 
The amplitude of the decay 7 — )■ a + /3 is 

Am{p, KA13) = {pK>^l3\U\jm) (5.3) 



^Helicity is defined as the projection of a particle's spin along its direction of motion, \ — S ■ P. 
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where C/ is an operator invariant under rotations and reflections. Using Eq. 5.1 this can 
rewritten as 

A^{p, Xa, X/s) = y^^M(A,, Xis)Dl^{^, ^, 4,) (5.4) 

where we have, for convenience, defined the matrix element M(Aq,,A^) = {jmXaXp\tJ\jm) 
which is independent of j and m. 

If the production mechanism of 7 is uncertain, the pure quantum state is unknown, and 
we can not calculate a scattering probability directly from this amplitude. Instead, we will 
use a density matrix to represent a more generalized quantum state. For any system of total 
angular momentum j and magnetic sub-state m, we can define a density matrix in the {j-,m) 
basis as 

p=^ \jm)pjn,m'{jm'\. (5.5) 

m,m' 

The density matrix is Hermitian and has a trace of 1, meaning that it can be characterized by 
4j(j+l) independent real numbers. If we force the x and z-axes to be in the production plane 
(the plane formed by the two incoming hadron momenta boosted into the 7 rest frame), we 
can further restrict the matrix with 

Pm,m' = (~lj P-m,-m'i I^-Oj 

which reduces the number of independent parameters to 4j. 

Using the density matrix we can now write the probability for 7 — t- a + /3 with p specified 
by Euler angles 99, ■!?, and ip in the 7 rest frame: 

^^4^ = ^^ Y^ Y^ Am{pXo,Xp)prarn'A^,{pXaXp) 
m,m' AqA^ 

= ^^Y Y.\^^^o.,Xp)\'D^^^{^,d,i^)Dl,^{^,^,i,)prnm'. (5.7) 

m,m' AqA^ 

For the specific case of 7 as a J/^, we have a spin-1 particle with m=0, and helicity conservation 
requires Aq / A/3. For decay leptons, A^ and X/s can be either | or — | so that the sum simplifies 
to 



d^qdn Aw 



<-^ 
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After requiring the x and z-axis to be in the production plane, the density matrix can be 
simphfied in terms of 4 real- valued parameters 



/ 



. \ 



/ 



Pii Re{pio) + i/m(pio) pi_i 

Re{pu)) - Um{pu)) 1 - 2pii -i?e(pio) + Um{pio) 

pi-i -~Re{pio) - Um(pio) pii 



Pii pio Pi- 
Poi poo po-i 

\y p_ii p_io p-1-1 y \^ 
and we can use the property c?^/^ = (— l)™"™- d^^, = d?_^_^i along with the contributing 
d-functions 






d\M 



1 + cos(i?) 



<<!„(*) = zn^ ,[_,(,) = 1;:^ 



to write out the full amplitude 



da 



d'^qdn 47r 



^2 2 



{pii(l + cos^ T?) + (1 - 2pii)(l - cos2 t?) 
+i?e(/9io) v2sin2'i? cos ip + pi-isin -d cos 2ip}. 



(5.9) 



It should be noted that the angular distribution has no dependence on Im{pio), implying that 
a measurement of each angular decay coefficient still allows for an undetermined phase in the 
density matrix. 

It is customary (76; 123) to parameterize the cross-section as 

1 2 

.2^ "2) {Wt{1 + cos^^) + Wl{1-cos'^^) 

+Wa V2sin2'd cos ip + WAAsin^'& cos2if} (5.10) 



da 3 

d^qdn ~ 4^ 



M( 



where Wt and Wl are called the 'transverse' and 'longitudinal' components respectively. The 
reason for the 'transverse' and 'longitudinal' labels is that Wt corresponds only to the pu and 
P-i-i elements of the density matrix (spin transverse to the i-axis), while Wl corresponds 
only to the /?oo element (spin along the z-axis). Plots of the decay distributions for purely 



longitudinal and purely transverse spin- alignment can be found in Fig. 5.4 Wa and Waa are 
called the 'single-spin flip' and 'double-spin flip' components. 

The angular distribution of the decay leptons is often rewritten as 



da 
d{cos'd)dip 



oc 1 + A^ cos^ -d + X-ffip sin 2'd cos 93 + X^sin'^'d cos 2ip 



(5.11) 
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wnere A^j = ^^^^^ - - ^_p^^ , a^^ = ^^+^^ - i-p^^ > ana A^ = ^^,+^1/^ - j^^- 
Experimentally, only the A^ parameters have typically been measured (21; 69; 4; 14), but 
this is clearly not adequate to fully describe the decay. Several recent experiments have also 
measured A^<^ and A<^ (3), and in Chapter [sl a new measurement of A^, A^<^, and A<^ will be 
presented from PHENIX data. 

The coefficients A^, A<^, and A^^ are clearly correlated. In particular, a rotation of the 
reference frame in the production plane by an angle 

5 = Jarctanf-^^) (5.12) 

2 V A(/3 — A^ / 

(90° when A^=A^) corresponds to a rotation into the frame in which A^;^ is zero, while A^ and 
A(^ describe the entire angular distribution (87). 

In hadronic collisions, intrinsic partonic transverse momenta allows asymmetries in the 
spin alignment to exist with respect to the production plane on an event-by-event basis. Such 



asymmetries modify Eq. 5.11 to include two additional coefficients. For a subprocess i. 



d{cos^)d^ "" ^ + ^^'^ '°'' "^ + ^^l ^''''^"' ""° ^ ^ ''^ 



ocl + Xf cos^ ^ + Xf^ sin(2i?) cos ip + X^^hiri^-d cos 2ip 

+Xi^'^ sin(2t?) sin(/7 + X^^'hin'^^sm2^. (5.13) 



For inclusive measurements, however, these additional coefficients must vanish (85), leaving 



the expression in Eq. 5.11 



5.4 Reference Frames 



The derivation of Section 5.3, leading to Eq. |5.11 gave angular distributions in terms of i? 



and (f with the z and x axes required to lie in the production plane. However, there are an 
infinite number of possible reference frames, each leading to an alternate definition of ?? and 
(p. Three of these reference frames are typically chosen for measurements of angular decay 
distributions: 

1. The Jacob- Wick Helicity Frame (111) 

The z-axis is chosen as the J/ip velocity direction in the lab frame. 
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Figure 5.4 Angular distributions of decay products from a spin-1 bo- 
son decaying into two spin-1/2 fermions for 'transverse' 
spin- alignment on the left (A^=l, A(^=A^^=0) and 'longitu- 
dinal' spin-alignment on the right (A^=-l, X^=X'd^=0)- The 
z-axis points towards the top of the page. 



2. The Gottfried-Jackson Frame (100) 

The i-axis is chosen as the 3-momentum of one beam boosted into J/ip rest frame. 

3. The Collins-Soper Frame (76) 

A = Pb Pa 

\Pb\ \Pa\ 

where Pa and Pb are the 3-momenta of each beam boosted into J/ip rest frame. 

After the z-axis is chosen, the y-axis is taken as normal to the production plane, y = z x Pfy, 
and the x-axis is x = y x z. A cartoon of the various frames and the angles associated with 



each frame can be found in Fig. 5.5 The Helicity frame has typically been used by collider 
experiments, while the Gottfried-Jackson and Collins-Soper frames have been used by fixed 
target experiments. The Collins-Soper frame was initially motivated by Drell-Yan production 
and has been widely used for measurements of that process. 

The orientation of the y-axis is especially important for the determination of X,^^. While 
the A^ and A;^ coefficients do not depend on the orientation of the y-axis, the sign of A^^ 
does. To make the situation even more complicated, the sign of A^,^ changes depending on the 



sign of the rapidity (see Fig. 5.4). In order to keep A^i^ from becoming zero when integrating 
over rapidity, the direction of the z and y-axes should be flipped when going from positive to 
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Figure 5.5 Illustration of z-axes for several reference frames (left) and for 
the angles -d and ip as defined in those frames (right). 



negative rapidity (86). 

One additional subtlety is that the Gottfried-Jackson frame is ambiguous in a collider 
environment, because either beam can be boosted into the J/ip rest frame to define the z-axis. 
This ambiguity is especially troublesome at forward rapidities where the detector acceptance 
is quite different depending on the beam used in the frame definition. To avoid ambiguity, 
we define two separate Gottfried-Jackson frames: the Gottfried-Jackson forward frame, using 
the beam circulating in the same direction as the J/^ momentum, and the Gottfried-Jackson 
backward frame, using the beam circulating in the opposite direction. 
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Figure 5.6 Orientation of the y-axis depending on the direction of the 
Lorentz boost. 
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Once a reference frame is defined and the various amplitudes are measured, it is possible to 
make calculations for other frames by rotating the density matrix (40). The fact that the real- 
space distribution of decays must be invariant under such rotations, along with the property 
that |A^| < 1 (by definition) lead to the inequalities 

|A^|<^(1 + A,,) , |A^^|<^(1-A^), (5.14) 

implying additional bounds |A(p| < 1 and |A^^| < 1 (84; 85; 86). These inequalities also allow 

for a family of frame-invariant quantities]^ one of which is 

~^A^^+3A^^ (5.15) 

1 — Xip 

which corresponds to the weighted sum of A^ over all n contributing subprocesses 

^ /« .« 



^3 + A«^'' 






4 = 1 -J + ^^ 

where /(*^ is the fraction of measured particles produced by process i. 

It should be noted that the entire formalism developed in Sections |5.3| and |5.4| is valid for 
any spin-1, m=0, particle decaying to two leptons, including both the J/ip and photon. A 
naive model for Drell-Yan production (the decay into two leptons of a virtual photon from 
the annihilation of a quark and anti-quark) would expect A^ = 1 and A^^=A(^=0, because 
real photons are transversely polarized. Virtual photons are likely to maintain much of that 
transverse polarization and transfer it to the decay leptonsjj However, non-zero total transverse 
momentum of the decay products can lead to A.,? ^ 1 (70), and non-zero transverse momentum 

of the colliding partons can lead to A^,^ / and A,^ / (72)n 

*The derivation of these quantities uses the property d\ ,„{■&) + dLi_m ~ ^|m|,i and assumes that a reference 
frame exists wherein only the pn and p_i_i density matrix elements are non-zero. 

^There is an unfortunate ambiguity in the use of the words 'transverse' and 'longitudinal' when describing 
polarization throughout the hterature. Typically, when the polarization of a gauge Boson is being discussed, 
as we are discussing the photon here, the labels correspond to the orientation of the field vectors (electric and 
magnetic in this case) with respect to the momentum of the gauge boson. For hadrons, however, the same labels 
usually correspond to the orientation of the hadron spin (which is perpendicular to the field vectors) with respect 
to the momentum of the hadron. When we discussed transversely polarized protons in Chapter |4] for example, 
the 'transverse' label meant that the proton spin, not the field vectors, was transverse to its momentum. To 
avoid this ambiguity as much as possible, we will not speak in this document of the ' J/V' polarization' as the 
spin-alignment is called elsewhere in the literature. 

®The A^ coefficient in Drell-Yan production is also sensitive to the Boer-Mulders effect discussed in Sec- 



tion 4.4 (52) 
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Lam and Tung have shown, however, that the spin-1/2 nature of the quarks leads to the 
constraint that 1 — A^ — 4A^ = 0, meaning that Wl = 2Waa (123; 124) F] While pion induced 
Drell-Yan experiments on deuterium and tungsten show clear violations of the Lam- Tung 
relation (102; 107), measurements of proton induced Drell-Yan on a deuterium target show 



no such violation (152). The expression in Eq. 5.16 is formally equivalent to the Lam- Tung 
relation when A=l, and any violation of the Lam- Tung relation could be described by a suitably 
modified frame-invariant expression (85). 

5.5 Production Mechanism 

Since the first observation of the J/^ meson, theorists have been attempting to model the 
mechanism for its production. Among the most prominent of these models have been the Color 
Evaporation Model, the Color Singlet Model, and the Color Octet Model. Each model makes 
distinct predictions for the production cross section and most also predict the angular decay 
coefficients of the produced charmonium states. 

If the mass of the charm quark is taken to be large compared to Aqcd, then a system of 
two hadronizing charm quarks is approximately non-relativisticj^ All three models attempt to 
factorize J/ip production into a relativistic part describing the production of the charm and 
anti-charm quark, dacc[n]+Xj and a non- relativistic part describing the bound state of the two 
quarks, Fcg[„](A): 

tZA^^^^Fcc-[n](A) (5.17) 

where the [n] denotes the quantum state of the cc pair and A the energy scale. 

5.5.1 Color Evaporation Model 

The Color Evaporation Model (CEM) (89) was the earliest attempt to calculate the char- 
monium cross section and has only been applied to hadronic collisions. The non-relativistic 



'^The Lam- Tung relation is analogous to the Callan-Gross relation in DIS, a consequence of the interaction 
betwee n a p hoton probe and half-integer spin quarks, resulting in the condition Fl — F2 — 2xFi, discussed in 
Section 



3.2 



*In the MS factorization scheme, Aqcd is approximately 206-231 MeV for 5 active flavors of quarks (54), 
whereas the charm quark rest mass energy is approximately 1.27 GeV (31). 
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part is assumed to be non-zero and constant between Arn^ and 4m^ and zero for all other 
energies, where rric and mo are the mass of the charm quark and D meson respectively. The 
model does not take the quantum state of the cc pair into account explicitly, but instead sums 
over all quantum states and multiplies by a factor of g, the statistical probability for the cc 
pair to be in the color-singlet state. The differential cross section is 

da{J/i; + X) = ^^ Y. n d^'^^-^^JW^^ (5.18) 

which has a single constant parameter i^ccfj/i/)! to be determined from a fit to data. 

The CEM is actually more of a parameterization than a model, but it has predictive power 
for the shape J/'ip cross-section to the extent that all dcr^cfnl+x can be calculated. There is 
no clear prediction for the J /ip spin-alignment from the CEM, but it has been suggested that 
multiple soft gluon exchanges destroy the spin-alignment of the cc pair (32). A comparison of 



the CEM to cross-section data shows impressive agreement (Fig. 5.7). 



5.5.2 Color-Singlet Model 

The Color-Singlet Model (CSM) (103; 50) is an attempt to explicitly take into account 
the quantum state of the cc pair. In this model, the cc pair emerging from the relativistic 
scattering diagram is assumed to be in the same quantum state as the produced J/V', and the 
non-relativistic amplitude is the real-space J/^ wavefunction evaluated at the origin: 



oo 



d(Tcc[iSi]+X 



da{J/ij + X)= dM -';;'"- V'j/^(r = 0). (5.19) 



dM 



The wavefunction ipj/^ is determined using hydrogen-like potential models for the hadronized 
cc pair (105), many of which reproduce the charmonium spectroscopy quite well, and the 
model has no free parameters. A comparison of a CSM calculation at next-to-leading order in 



as (NLO) with data from RHIC is shown for the differential cross section in Fig. 5.8 and the 



angular decay coefficient A^ in Fig. 5.9 (125). While the measurement of Xg does not constrain 
the model very well, the CSM clearly under-predicts the J/ip cross section at both mid and 
forward rapidities. NLO diagrams increase the cross-section relative to leading order (LO), 
but the increase is not large enough to remove the discrepancy. 
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Figure 5.7 Comparison of Color Evaporation Model predictions to mea- 
sured cross-sections for production of direct J/^p (top left), 
prompt J/ip from decays of V'' (top right), and prompt J/tp from 
decays of Xc (bottom) at CDF as a function of px- The dashed 
and sohd curves use the MRST98 HO (higher order) (129; 128) 
and GRV98 HO (98) PDF sets with mc=1.2 and 1.3 GeV/c^ 
respectively (from (60)). 



Similar conclusions have been drawn about the CSM from comparisons of calculations at 
NLO with data from HERA and Zeus (45; 68). Recent calculations at NNLO, however, seem to 
accurately reproduce the measured T cross-section at CDF (44) . The slow convergence of the 
perturbation theory seems to imply that the model does not properly order the contributing 
diagrams, and it is quite possible that important effects are either not yet included or are not 
included at the correct order (126; 134; 106; 43). 



5.5.3 Color-Octet Model 

The Color-Octet Model (COM) (56) attempts to formalize the factorization of relativistic 
and non-relativistic effects. The model uses a generic expansion 



da{J/^ + X) = Y, 



dM- ,,, 
dM 



-{0, 



(5.20) 
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Figure 5.8 Comparison of the CSM at NLO in as with the inclusive J/^ 
cross-section from PHENIX and STAR at mid rapidity (left), 
and with the PHENIX inclusive J/ip cross-section at forward 
rapidity (right) (from (125)). Both plots assume that 59ibl0% 
of the J/ip mesons are directly produced, and in both cases the 
model under-predicts the data. 



with parameters (Or^i ); non-relativistic matrix elements associated with the amplitude for 
producing a J/ip from a cc pair in state [n]. Techniques developed in non-relativistic QCD (66) 
are then applied to determine the size of the {0,J^ ) parameters in powers of v, the relative 
velocity between the c and c, where f'^~0.3c^ for the J/ip (121). The model is thus a double 
expansion, about v'^ and as- The COM predicts that the leading color-singlet diagram is of 
order a^(2?n,c)^/p|i, while the color-octet is of order a^{2mc)'^v^/pj^. For J/tp mesons from 
fragmentation, the leading color-singlet diagram is of order a^/p'^, and the color-octet is of 
order a^v^/p"^ (122). These diagrams are shown schematically in Fig. 5.10 and 5.11 At 
large pT, color-octet fragmentation becomes especially important, leading to the predicition 
of transverse spin-alignment at large pT in the Helicity reference frame for the COM (because 
the majority of the J/ip spin should come from a single transversely polarized gluon) (53; 59). 
The COM takes its name from the prediction that color-octet diagrams dominate for J/ip 
production at large pT for hadron colliders like the Tevatron, but it also predicts that color- 
singlet diagrams dominate for J/ip mesons with smaller pT from e^e^ collisions at B-factories 
like KEKB. A comparison of a COM calculation at next-to-leading order in ag (NLO) with data 



from CDF is shown for the cross-section in Fig. 5.12 and for the A^ angular decay coefficient 
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Figure 5.9 Comparison of the CSM at NLO in as with the a spin-aHgn- 
ment coefficient (called 'A^' in this document) for inclusive J/^ 
in the Helicity frame from PHENIX at mid rapidity (left) and 
forward rapidity (right) (from (125)). The data does not con- 
strain the model very well. 
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Figure 5.10 Feynman diagrams for color-singlet charmonium production 
at leading-order (left) and from color-singlet fragmentation 
(right). 



in Fig. 5.13 (99). The model agrees very well with the cross-section but drastically disagrees 
with A^. 



5.6 Discussion 



None of the J/tp production mechanism models discussed here adequately describe the 
data. Charmonium cross sections from hadronic collisions are well matched by the CEM, but 
the model does not make predictions for DIS or e~^e~ collisions and has very little predictive 
power for angular decay coefficients. The CSM makes predictions for the A^ angular decay 
coefficient which agree with data, but the model severely underpredicts measured J/ip cross 



46 



TfTOWr 



viUUUU£ 





Figure 5.11 Feynnian diagrams for color-octet charnionium production for 
leading-order production in the t-ctiannel (left) and from col- 
or-octet fragmentation (right). 
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Figure 5.12 Comparison of the COM at NLO in as with the cross-sec- 
tion for directly produced j/ijj mesons from CDF as well as 
predictions for the cross-section at the LHC (from (99)). 



sections. The COM predictions, meanwhile, match measured cross sections but disagree in 
both sign and magnitude with the A^ angular decay coefficient. 

Unfortunately, there is also some tension between separate measurements of A^. CDF, the 
leading experiment reporting A^ at large pT, found A^ consistent with the COM from Run-I 
data (21) but inconsistent with the model from Run-II (4). While the change of conclusion 
has been attributed to the identification of systematic effects which were not accounted for in 
the analysis of the Run-I data, further measurements are necessary to confirm disagreement 
with the COM. 

Previous measurements of A^ at mid-rapidity from p + p collisions at y^=200 GeV from 
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Figure 5.13 Comparison of the COM at NLO in as with the a spin-ahgn- 
ment coefficient (called 'A^' in this document) for directly pro- 
duced J ftp mesons from CDF (from (99)). 



PHENIX (14) are limited to small px, where predictions have been made from both the 
CSM (125) and COM (71), but the COM predictions are questionable at small pT with- 
out further theoretical effort. Furthermore, the data are not able to distinguish between the 
predicted A^ from the COM and CSM in the region where the calculations have been made. 

A change in center of mass energy from ^/s=200 GeV to ^/s=500 GeV increases the J/ip 
cross section (92). An increase in y/s, along with the larger coverage of the PHENIX muon 
spectrometers relative to the central spectrometers (2 units of rapidity as opposed to 0.7), mean 
that a measurement at forward rapidity from y/s=500 GeV p + p collisions at PHENIX allows 
for significantly higher statistics. Measurements with larger statistics at forward rapidity, 
which will be presented in Chapter [8j bring the pT coverage of PHENIX to a region where the 
COM is less questionable. 

Furthermore, the full azimuthal coverage of the PHENIX muon spectrometers allows for 
measurements of the X^p and A^<^ coefficients in addition to A^. Measuring all three coefficients 
in several reference frames gives a much deeper understanding of J/^p production, providing 
several reference points for comparisons to production mechanism models. Recent measure- 
ments from p + N collisions at -^8=41. 6 GeV from HERA-B (3) suggest that the Collins-Soper 
z-axis is closer to the natural axis for J /ifj spin alignment than the i-axis of the Helicity frame, 
and theoretical attempts to determine the true natural axis have already begun (58). Mea- 
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surements of all angular decay coefficients from PHENIX in several reference frames, as I will 
present in Chapter [8j can either confirm or contradict observations from HERA-B and provide 
input to theory, leading to a more comprehensive understanding of J/^ production. 
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CHAPTER 6. RHIC and the PHENIX Experiment 



6.1 RHIC 

The Relativistic Heavy Ion Collider (RHIC) is a circular accelerator with 0.6 km radius in 
Upton, NY capable of colliding heavy ions (Au, Cu, Si, etc.) at beam energies up to 100 GeV 
and polarized protons at beam energies up to 250 GeV. 

For polarized proton collisions, an optically pumped polarized H~ source provides protons 
to a Linear Accelerator (LINAC), which accelerates them to 200 MeV and injects into the 
Alternating Gradient Synchrotron (AGS). The AGS further accelerates the protons to 1.5 GeV 
before injecting them into the RHIC accelerator ring, where they are brought to full energy 
and eventually into collision. 

Beam polarizations are measured by two independent polarimeters, a fast carbon-target 
polarimeter (133) for relative polarization measurements and a hydrogen-jet polarimeter (135; 
83) for an absolute measurement. The carbon-target polarimeter measures the polarization 
for each store of beams in the RHIC accelerator, and the hydrogen-jet polarimeter is used over 
a longer time-period to obtain a calibration for the carbon-target measurements. The stable 
polarization direction in the RHIC ring is transverse to the direction of the protons' motion, 
and the polarimeters take advantage of the transverse single-spin asymmetries (SSA) discussed 
in Section [4] in order to measure the beam polarization. 

For the hydrogen-jet polarimeter, a polarized hydrogen jet is brought into collision with each 



of the proton beams, and recoil protons are measured by silicon detectors (see Fig. 6.1). Data 
from the detectors is used to calculate Ajy ^target, the transverse SSA with an unpolarized beam 
and polarized target, and Ajy h^^arm the transverse SSA with a polarized beam and unpolarized 
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target jj The beam polarization is then 



Vbe 
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target 



where the target polarization Vtarget is independently measured 




(6.1) 



Figure 6.1 Setup of the RHIC hydrogen-jet polarimeter. Silicon detectors 
on the right and left measure the transverse single spin asym- 
metry of the recoil protons. 



In the carbon-target polarimeter, an unpolarized ultra-thin carbon target is brought into 
collision with the polarized proton beams. The transverse single spin asymmetry j47v,6eam is 
then determined in each store by measuring the spin-dependent azimuthal distribution of scat- 



tered ions with silicon detectors surrounding the target (Fig. 6.2). Because the transverse SSA 
scales with p, only a single absolute measurement of "P (from the hydrogen-jet polarimeter) is 
necessary to calibrate the carbon-target polarimeter and determine polarizations for all stores. 



6.2 PHENIX 



The Pioneering High Energy Nuclear Interaction eXperiment (PHENIX) is a large, multi- 
purpose experiment on the RHIC accelerator ring. It was designed with fine granularity, 
leading to good position and momentum resolution at the sacrifice of acceptance to produced 



^In reality both the beam and the target are always polarized, and statistics are combined in such a way 
that the asymmetry is only sen sitive to one polarized particle. The procedure and consequences for doing this 
will be discussed in Section 7.3 
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Figure 6.2 Setup of the RHIC carbon-target polarimeter. Silicon detec- 
tors surround the target (represented by the cross-hairs in the 
center of the image) . The proton beam momentum direction is 
through the page. 



particles (16). The experiment consists of four large spectrometers: two for tracking and 
calorimetry at central rapidity and two for muon tracking and identification at forward and 
backward rapidities. The central detectors cover |?7| <0.35 and A(j) = 2 x |, while the muon 
detectors cover approximately 1.2< |r/| <2.2 and A0 = 2tt. For global event characterization, 
there are beam-beam counters (BBC) to measure event time and position. Additionally, a 
local polarimeter is used to monitor the proton polarization direction. A schematic view of the 



various detector subsystems can be seen in Fig. 6.3 
2008 A PHENIX Detector 
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Figure 6.3 Drawing of the PHENIX Detectors in beam view (left) and side 
view (right). 
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6.2.1 Detectors for Event Characterization 

In order to make a useful measurement with a high energy collider, it is typically necessary 
to determine both the time and position of each collision. A determination of the time of 
the collision is necessary in order to tag the crossing in which the collision occurred, and 
the event vertex is necessary for tracking. A measurement of the collision time is especially 
important for spin-dependent measurements, where the collision must be restricted to a single 
bunch in the accelerator in order to identify the polarization direction of the colliding protons. 
PHENIX uses beam-beam counters (BBC) to measure the time and longitudinal position of 
each cohision (28) j^ 

The BBC consists of two identical counters positioned 144 cm on either side of the center 
of the experiment. Each counter consists of 64 photo-multiplier tubes (PMTs) mounted on 
a 3 cm thick quartz radiator in a cylinder of outer radius 30 cm and inner radius 10 cm 
(corresponding to approximately 3.0 < \rj\ < 3.9). Raw signals from the PMTs are sent 
through time-to- voltage converters to flash analog-to-digital converters, which pass digitized 
timing and pulse-height information to a local level-1 trigger (LLl). The LLl makes a decision 
regarding whether or not the raw signal information should be sent to an event builder and 



recorded (triggering will be discussed in more detail in Section 6.2.5). The RHIC beam clock 
operates at approximately 9.6 MHz so that a crossing occurs every 104 ns. Each BBC has an 
RMS timing resolution of 54ib4 ps, which is clearly fine enough resolution to determine the 
crossing in which a collision has occurred. 

6.2.2 Local Polar imeter 

While RHIC is responsible for measuring beam polarization using the polarimeters dis- 
cussed in Section |6.1[ PHENIX is responsible for monitoring the direction of that polarization 
near the collision vertex. This monitoring using measurements of neutrons from a zero-degree 
calorimeter (ZDC) (17) and shower maximum detector (SMD), collectively referred to as the 
local polarimeter. 



■^The transverse size of the colliding beams is less than 100 ^m and has limited effects on reconstructed 
particle momenta. 
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The ZDC consists of three layers of hadronic calorimeter located approximately 18 m from 
the center of the interaction region. Each layer is composed of a Cu-W alloy absorber and a 
Hamamatsu R329-2 photo-multiplier, which collects Cerenkov light from optical fibers behind 
the absorber jj The entire layer corresponds to 1.7 interaction lengths of material, and the 
detector has an energy resolution of 21% for a 100 GeV neutron. An additional 3.3 mm layer 
of scintillator is placed in front of the ZDC for the identification and veto of charged particles. 
Neutrons are typically selected by requiring an energy deposit between 20 and 120 GeV along 
the beam direction in the ZDC along with less than one minimum ionizing particle in the veto 
scintillator. 

Interleaved between the first and second layer of the ZDC are position sensitive SMDs 
composed of 7 15 mm wide horizontal strips and 8 20 mm wide vertical strips. The SMDs 
are tilted at an angle of 45° away from the interaction region with a position resolution of 
approximately 1 cm and an active area covering ~0.3-1.4 mrad from the beam. 

To determine a polarization direction, the transverse SSA of forward neutrons is measured 
as a function of azimuthal angle from the nominal spin direction and fit with 

A{4>) = An sm{^-M (6.2) 

to monitor a phase (po from the nominal spin orientation. 

6.2.3 Central Arm Detectors 

The PHENIX central arm spectrometers are composed of separate subsystems for tracking, 
particle identification (PID), and calorimetry: 

1. Calorimetry is done using Pb-scintillator sampling calorimeters and a Pb-glass Cerenkov 
calorimeter. (39). 

2. PID for electrons is handled by a ring-imaging Cerenkov detector (24). 

3. Tracking uses drift chambers and pad chambers (15). 



''Since neutrons are not charged, they cannot directly produce Cerenkov radiation. Instead, the radiation is 
produced by fast electrons emitted as fission products from neutron interactions in the calorimeter. 
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All contributing subsystems are used in coordination to make a meaningful measurement (132). 

6.2.3.1 Central Arm Calorimetry 

Two subsystems located at the outer radius of the PHENIX central arms make up the 
electromagnetic calorimeter (EMCAL). Approximately 75% of the azimuthal acceptance of 
the detector uses a Pb-scintillator sampling calorimeter (PbSc), and the other 25% uses a Pb- 
glass Cerenkov calorimeter (PbGl). The PbGl has better granularity and energy resolution, 



while the PbSc has better timing resolution (Table 6.1), making the two subsystems quite 
compliment ar y. 



A diagrammatic drawing of a PbSc module is shown in Fig. 6.4 Each module consists of 
four towers composed of 66 5.535 cmx 5.535 cm sampling cells of alternating Pb and scintillatoi|j 
with edges plated in Al. The Pb acts as a passive absorber, while the scintillator samples 
incoming energy. Light is collected by 36 fibers in each cell and read out by photo-tubes in the 
back of the tower. One 'super-module' is made up of 32 modules mounted on a single stainless 
steel skin. Eighteen of these super- modules are combined to make one 2x4 m^ sector with its 
own steel support frame. 

The PbGl section of the central calorimeter is composed of 192 super-modules, each with 
24 separate Pb- glass modules. The modules have a high index of refraction so that high 
energy charged particles moving through them radiate Gerenkov light, and each module has 



a dedicated photo-multiplier tube for readout. Fig. 6.5 shows an exploded view of a PbGl 
super-module along with the LED board consisting of dedicated LEDs for calibrating each 
module. 

Table 6.1 Energy, spatial, and timing resolution of the PHENIX PbGl and 
PbSc calorimeters for electrons. 

Detector aE{E)/E cFx{E) (mm) at{E) (ns) 

8.1% 1 t^r:^ 5.7 n n^ ^ 0.03 



PbSc 2.1% e^M^ 1.55 // 0.06 u. p,r,v nm 
^E/GeV ^E/GeV g/GeV-O.Ol 

PbGl 0.8% e -^MZ 0.2 e , ^-^ i 0.075 e^J^li^ 

yjE/GeV ^E/GeV ^500-E/GeV 



''The scintillator used in the sampling cells is p-bis[2-(5-Phenyloxazolyl)]-benzene (POPOP) with a fluorescent 
additive p-Terphenyl (PT) 
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Figure 6.4 Diagrammatic view of a PHENIX Pb-scintillator calorimeter 
module. 



6.2.3.2 Central Arm Particle ID 

Located between the calorimeters and tracking detectors is the ring-imaging Cerenkov 
(RICH) detector, used for discriminating between electrons and pions below the pion Cerenkov 
threshold of 4 GeV/c. The majority of the RICH volume consists of 40 m^ of ethane held at 
0.5" of water above ambient between two glass plates. Electrons moving through the ethane 
radiate a cone of Cerenkov light, which is focused into ~14.5 cm rings by spherical mirrors 
onto an array of photo-multiplier tubes. Pions with momenta below 4 GeV/c do not produce 
Cerenkov light so that no ring is formed, and they can easily be rejected. Imperfections in the 
focusing mirrors lead to a 2.5 mm inaccuracy in the ring position. 



6.2.3.3 Central Arm Tracking 

The Drift Chamber (DC) is the innermost detector of the central arm spectrometer and 
provides the primary momentum measurement and vector for tracking. It is outside of the 
primary magnetic field, but a residual field of approximately 0.6 kCa accounts for ~1° of track 
bending within the chamber. The detector consists of 2m long cylindrically shaped detectors 
located 2-2.4 m from the beam axis. Each chamber consists of 20 sectors, made up of 6 wire 
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Figure 6.5 Exploded view of a PHENIX Pb- glass calorimeter super- mod- 
ule. 



modules each. There are a total of 6,500 anode wires, making up 13,000 readout channels, and 
leading to 165 /um single wire resolution and approximately 2mm spatial resolution on each 
track. 

Outside of the DC are 2 layers of non-projective pad chamber (PC), one directly inside of 
the RICH and one inside of the EMCAL. The West arm has an additional layer of PC directly 
outside of the RICH. The inner layer of PC provides a z-coordinate at the exit of the DC and 
assists in projecting tracks through the RICH, while the larger outer layer provides matching 
to clusters in the EMCAL. Each PC is segmented into pixels 8.2 mm in z by 8.4 mm in i-(j). In 
order to consolidate readout channels and minimize noise, 9 pixels are ganged into one 'pad' 



and hits in 1 cell (3 pads) are required to register a single hit (see Fig. 6.6). 



6.2.4 Muon Arm Detectors 

At forward rapidities, spectrometry is devoted mainly to muon detection by two major 
subsystems, the Muon Tracker (;uTr), used for tracking, and the Muon Identifier (juld) for 
muon identification (25). The trajectory of a muon entering the detector is bent in azimuth by 
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Figure 6.6 PHENIX Pad Chamber pad and cell configurations. On the left 
is one pad, composed of 9 pixels, and on the right is a number 
of interleaved pads. The outlined boxes correspond to a hit in 
one cell (3 pads). 



a radial magnetic field of several kGa, shown schematically in Fig. |6.8| (magnets are described 
much more fully in (41)), so that their momenta can be measured by the finely segmented 
/xTr. The majority of hadrons are absorbed by the thick steel plates of the /ild so that a muon 
trigger can be formed in the furthest gaps of the detector from the collision vertex. 

On each side of the collision vertex are 3 /iTr stations with 3 chambers in the closest two 
stations and two in the furthest station from the vertex. Each chamber is constructed in octants 
with the 5 mm wide cathode strips positioned at 11.25° from radial. The chambers are filled 
with a gas mixture of 50%Ar + 30%CO2 + 20%CF4 and held at a voltage of approximately 
1850 V. As an ionizing particle moves through the gas, ~80 fC of charge is deposited on the 
cathode strips. 

Behind the ;uTr on the north side of the experiment is a 30 cm thick steel backplate for 
the magnet, and a 20 cm thick backplate is behind the /uTr in the south. The backplates are 
followed by 2 10 cm and 2 20 cm steel absorbers, creating a total of five gaps. The gaps are 
instrumented with four planes of 8.4 cm wide larocci tubes each, two oriented vertically and two 
horizontally. Signals between the two sets of tubes in the same orientation are OR'ed together 
creating a signal with higher efficiency and faster signal timing than would be produced by a 
single tube. A muon must have an energy of 2.7 GeV in order to reach the final gap of the 
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Figure 6.7 Cut-away view of the PHENIX muon detectors. The three sta- 
tions of the /uTr are located in the yehow structure while the 
;uld planes are drawn in green. 



/xld, and the probability for a 4 GeV/c pion to reach the final gap is 3% or less. 



6.2.5 Triggering 

The PHENIX data acquisition system is limited to ~5 kHz of bandwidth for data through- 
put from all detectors, meaning that signals from subsystems cannot be constantly recorded. 
Instead, data must be 'triggered' using some minimal criteria for the detector signals. In 
PHENIX, this is done with several local level-1 (LLl) trigger systems combined into a single 
global level-1 (GLl) trigger. LLl triggers are typically specific to a single subsystem, and the 
GLl trigger combines trigger decisions from the separate LLl triggers. 

In order to avoid writing useless data to disk, it is essential to determine that a collision 
has taken place. The LLl system responsible for determining a valid collision time and vertex 
is the BBC LLl trigger, which takes raw input from the BBC and converts it into timing 
signals. Dead or noisy channels are masked in both arms, and the number of hits in each 
arm are determined. The trigger can then make a decision based either on the multiplicity in 
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Figure 6.8 Drawing of the magnetic field orientation in the PHENIX de- 
tector during the 2009 RHIC Run. 



each BBC (at least one hit in each detector signifies a possible collision) or select events more 
strictly by determining whether the collision vertex occurred within some distance about the 
center of the experiment. 

Triggering on electrons in the PHENIX central spectrometers is done using the EMCAL- 
RICH trigger circuit (ERT). The electron trigger of the ERT takes raw signals from the central 
arm calorimeters and sums the energy in either 2x2 or 4x4 squares of towers. If the energy 
summed within that region is above a specified threshold and there is corresponding activity in 
the RICH, the event is assumed to include a candidate electron and is triggeredjj In order to 
increase the rejection (number of events vetoed), the ERT is typically AND'ed with the BBC 
LLl trigger (e.g. both triggers must accept the event in order for it to be recorded). 

For muons, the /xld is currently used as the primary triggering system, although upgrades 
are underway to implement a momentum-sensitive trigger using the fiTr. The /ild LLl trigger 
uses raw signals from logical combinations of tubes containing all tubes within a legitimate 
track projection in the detector (tracks should be oriented such that they are coming from the 
collision vertex). Requirements are then made on the gaps hit within these projections. The 
trigger for 3/tp mesons, for instance, requires two separate logical units to have hits in 3 out of 



^Note that a photon trigger can easily be made by requiring energy above a given threshold in a number of 
towers with no requirement on the RICH. 
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5 gaps, one in either the first or second gap and another in either the fourth or fifth (furthest 
from the collision vertex). A single event meeting these requirements is called a 'deep' trigger, 



and a schematic for the logic is shown in Fig. 6.9 Occasionally, a similar 'shallow' trigger 



(Fig. 6.10), which only makes requirements on the first three gaps, is used in coincidence with 



a deep trigger to select J /ip mesons. Due to the poor timing resolution of the /ild, the ^Id 
LLl trigger is only used in coincidence with a BBC LLl trigger. 
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Figure 6.9 Logic for a 'deep' trigger in the /ild LLl. The A, B, and C 
represent separate groups of tubes contained in the logical unit 
for each gap. 



6.2.6 J/^ mesons at PHENIX 

The spectroscopy for electrons at mid-rapidity and muons at forward rapidity make PHENIX 
an ideal experiment for measuring J /ip meson decays to leptons over a large kinematic range. 
Extensive measurements of J/ij: mesons have been made at PHENIX from various collision 
systems, including: 

• Cross-section from -^5=200 GeV p+p collisions at mid and forward rapidity (8) 



Angular decay coefficients from y^=200 GeV p+p collisions at mid-rapidity (14) 
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Figure 6.10 Logic for a 'shaHow' trigger in the fjld LLl. The A, B, and 
C represent separate groups of tubes contained in the fogical 
unit for each gap. 



• Cross-section from -y/s/viv=200 GeV Au+Au colhsions at mid and forward rapidity (9) 

• Cross-section from ultra-peripheral ,/snn=200 GeV Au+Au collisions at mid-rapidity (20) 

• Cross-section from ^SN]\f=200 GeV d+Au collisions at mid and forward rapidity (19) 

• Cross-section from -y/s/v7v=200 GeV Cu+Cu collisions at mid and forward rapidity (10) 

In Chapter^ I will present a measurement of the J/ip transverse SSA from ^/s=200 GeV p+p 
collisions at both mid and forward rapidities, for the first time taking advantage of the polarized 
protons available to PHENIX in a J/ip measurement, and using nearly all of the subsystems 
discussed in this chapter. The measurement of Chapter [8] will present the first measurement 
of the J/^ angular decay coefficients using the muon spectrometers from ■y/s=500 GeV p+p 
collisions. 



62 



CHAPTER 7. J/ip Transverse Single Spin Asymmetry 



7.1 Motivation 

It was recently proposed (151) that the transverse single spin asymmetry (SSA) for direct 
J/ip production is sensitive to both the Sivers effect (Section |4.2[) and the J/ip production 



mechanism (Section 5.5). The direct production of J/'0 mesons is dominated by processes 



involving only gluons in the initial state, and as I stated in Section 4.3, gluons do not carry 
transversity. This implies that neither transversity nor the Collins or Boer-Mulders effects can 
contribute, and only the Sivers effect can be responsible for any transverse SSA. 

In the TMD approach used in (151), the transverse SSA is only non-zero if there is a non- 
zero gluon Sivers function and J/V' mesons are produced from p+p collisions in the color-singlet 
configuration, where only initial state radiation contributes. The transverse SSA cancels to 
all orders between the diagrams for initial and final state radiation for J /ip mesons produced 



in the color-octet configuration (151) (see Fig. 7.1). It is important to remember, however, 
that factorization is broken for the TMD approach to the Sivers effect in p+p collisions, es- 
pecially in the case of J /ij: production at RHIC, where the transverse momentum of the J /tp 
is approximately equal to its mass. The relationship between the transverse SSA and produc- 
tion mechanism is not as simple in the collinear approach, where there is partial but not full 
cancellation of transverse SSA generating terms for the color-octet configuration (118). 

Transversely polarized protons were brought into collision for approximately 4 weeks at a 
center of mass collision energy of y^=200 GeV during both the 2006 and 2008 RHIC runs, 
and both of these data samples are used to determine the transverse SSA for inclusive J/'0 
productionjj 

^ Since the measurement is made for inclusive J /ip mesons and not only for those which are directly produced, 
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(a) 



(&) 



(c) 



Figure 7.1 Diagrams which potentially produce an transverse SSA for J/^ 
meson production. Only initial state interactions contribute to 
color-singlet production (a). For color-octet production, inter- 
actions in the initial state (b) and final state (c) cancel in the 
TMD approach, causing the transverse SSA to vanish. (Fig. 
from (151)) 



In order to utilize the entire PHENIX acceptance, separate analyses are done for J/iJj^-e^e~ 
in the central spectrometers and J/tp^fi'^fi^ in the muon spectrometers. The analysis for the 
central spectrometer is done using data from the 2006 run and for the muon spectrometer 
using data from the 2006 and 2008 runs. Each analysis requires different considerations due 
to differences in the detectors and the number of J/ip mesons measured by each detector. 
I was primarily responsible for determining the asymmetry from the central spectrometers 
but also provided assistance in the final stages of the analysis in the muon spectrometers. 
For the benefit of the reader and for a more comprehensive view of the result, I will present 
measurements from both detectors. The result has recently been posted to the arXiv and 
submitted for publication (13). 



7.2 J/i) Selection 



To ensure that tracks used in the analysis are from electron and muon pairs from J/ip 
meson decays, cuts are applied to single tracks and to the pair. 



the theoretical guidance given above will be helpful, but not complete. A more complete theoretical calculation 
will need to take into account J/^p mesons from Xc decays. 
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For electron candidates, a coincidence is required between the BBC and the ERT triggersjj 
Electron tracks are required to have total momentum >0.5 GeV/c, to avoid background from 
electron conversions within the detector. The energy measured in the EMCAL divided by 
momentum measured in the DC is required to be within 4 standard deviations (cr) of 1 to 
ensure that the track is made by an electron, and the position matching between the track in 
the DC and the energy cluster in the EMCAL is required to be < 4cr in both beam direction 
and azimuth. Only events with a collision vertex measured by the BBC of ±30 cm about the 
center of the detector are considered, because the tracking of the central spectrometer does 
not work very well for collision vertices outside of that range. 

Muon candidates from the 2006 run are required to have a coincidence between the BBC, 



a fjID deep, and a filD shallow trigger (discussed in Section 6.2.5). For the 2008 run, a 
coincidence is required between the BBC and two /ilD deep triggers. The J /tp^ii^ ^^ analyses 
required the collision vertex to be within 35 cm of the center of the experiment along the 
beam direction]^ Each track is required to have momentum along the beam direction 1.4 < 
Pz{Ge\ /c) < 20 to avoid the trigger turn-on region, and the distance between the projection of 
the track in the ;uld and the track in the /iTr to the first ^Id gap is required to be <25(30) cm 
for the spectrometer on the North (South) side of the experiment. The angle between the 
projected tracks is required to be <10° in both the North and South detector. A fit of track 
pairs to the BBC vertex is required to have a x^ /ndf < 5 with 4 degrees of freedom. 

A plot of the invariant mass for electron pairs passing all cuts is shown in Fig. |7.2| and 
for muons in Fig. 7.3 The resolution of the J /ijj mass peak is found to be 61 it 2MeV/c^ in 



the central spectrometer and roughly 143 it 3 (154 it 3) MeV/c'^ in the North (South) muon 
spectrometer. 

In the central spectrometer, the contribution of pairs from random uncorrelated track 
combinations are estimated to be twice the geometric mean of the positive and negatively 
charged pairs, 2y/N^+f,+ N^-f,-, and subtracted from the spectrum of oppositely charged pairs. 



6.2.5 



^Descriptions of these triggers are given in Section i 

^A cut on the vertex is made at 30 cm by the BBC trigger, but the trigger does not perform a fuh vertex 
reconstruction, and the looser cut at 35 cm allows the inclusion of events with errors in the vertex determined 
by the trigger. 
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Figure 7.2 Invariant mass spectrum for electron candidate pairs and un- 
correlated track pairs in the central spectrometer. 



leaving only those pairs from collision-related processes. The fraction of remaining pairs coming 
from sources other than J/^ mesons is determined in a separate analysis using a full detector 



simulation of all physical processes contributing to the continuum (Fig. 7.4). Background 
fractions, defined as the ratio of correlated pairs from other physical sources to pairs from J/ip 
meson decays in a mass range from 2.7 < M(GeV/c) < 3.4 were estimated with the same 



simulation (Table 7.1). 



In the muon spectrometers, the number of background pairs is determined by fitting the 
invariant mass distribution of oppositely charged muon pairs with a function 



dN 
dM 



+ 



(ao + aiM + a2M^ + a^M^) 

e 2^ H ^e 



(A/-M^,,)^ 



27Ty/a 27rVcr' 

where the first and second Gaussian distributions approximate the shape of the J/ip and 
■0' peaks respectively, and the third order polynomial parameterizes the shape in invari- 
ant mass of both the uncorrelated pairs and the background from Drell-Yan production, 
p+p— T-DDX— T-yU^/i^X (open charm), and p+p^HJiX^fi^ n^ X (open bottom). Background 
fractions from the muon spectrometers, defined as the total number of background pairs di- 
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Figure 7.3 Example invariant mass spectrum from the 2006 dataset in the 
North muon spectrometer (p^ >0) for (a) oppositely-charged 
muon candidate pairs and (b) charged candidate pairs with the 
same sign. The solid line in the figure is the sum of Gaussian dis- 
tributions for the J/V' (dashed curve), ^' (dotted curve), and a 
third-order polynomial background (dotted-dashed curve) . The 
shaded regions show the area used for determining the back- 
ground Ajsj. 



vided by the number of signal pairs in a mass range 2a about the J/'0 Gaussian were estimated 



from the fit (Table 7.1). The background fractions are found to be larger for the 2006 dataset 



because the trigger requirement is not as restrictive, allowing for the reconstruction of more 
random track combinations. 

For the analysis in the central spectrometer, the number of e^e^ candidate pairs with 
invariant masses 2.7 < M(GeV/c^) < 3.4 divided by the number of counts in the BBCP] are 



"^The number of counts in the BBC is proportional to the provided luminosity, meaning that the number of 
J/^ mesons divided by the number of BBC counts should not change over time as long as the active area of 
the detector and the contributions from backgrounds do not change. 
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Figure 7.4 Invariant mass spectrum of electron pairs showing the contri- 
butions from ah physical sources. 



found to be constant over the course of the running period, implying that both the backgrounds 



and acceptance are approximately constant over time. (Fig. 7.5). Kinematic distributions for 



the accepted pairs are show in Fig. 7.6 



7.3 Asymmetry Formulae 



The general expression for a transverse single spin asymmetry (SSA) is given by Eq. 4.1 



Since both proton beams are polarized at RHIC, however, we never measure a single spin 
asymmetry directly using this expression but instead average over the polarization of one 
beam, defining the single polarized cross section as 



at = a^ + aU 



a^ = a'^ + a'^^ 



(7.1) 



where the f and ^ represent the polarization directions of the two beams. In the following 
discussions, t will be used for protons polarized radially outward from the RHIC ring and 4- 
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Table 7.1 Background fractions for J/ip mesons in the central and muon 
spectrometers at several transverse momenta. The background 
fraction for the central spectrometer includes only physical pro- 
cesses, while the background fraction for the muon spectrometers 
also includes random track combinations. 



PT (GeV/c) 


dataset 


detector 


background fraction 


0-6 


2006 


South Muon 


21.7±0.6% 




2006 


North Muon 


19.1±0.4% 




2008 


South Muon 


16.4±0.2% 




2008 


North Muon 


14.2±0.2% 




2006 


Central 


6.6±0.4% 


0-1.4 


2006 


South Muon 


23.2±0.7% 




2006 


North Muon 


22.0±0.7% 




2008 


South Muon 


16.1±0.3% 




2008 


North Muon 


15.5±0.3% 




2006 


Central 


5.6±0.5% 


1.4-6 


2006 


South Muon 


20.1±0.8% 




2006 


North Muon 


14.1±0.5% 




2008 


South Muon 


15.6±0.4% 




2008 


North Muon 


10.5±0.2% 




2006 


Central 


7.8±0.7% 



for those polarized radially inward for the analysis of 2006 data. For the analysis of 2008 data, 
they denote up and down. Using Eq |7.1| we can perform statistically independent analyses for 
each polarized beam and combine the results to get a smaller statistical uncertainty in the final 
result. 



Combining Eq. 4.1 and Eq. |7.1| gives 



A 



N 






(7.2) 



where A^^t^ N^-^^ N^^ ,ai\.(i N^^ are the experimental yields in each spin configuration, and 
TZi = C^ / C^^, 7^2 = C^ / C^"^ , and TZi = C^ / C^^ are the ratios of provided luminosities. By 
using a single polarization "P, we have implicitly assumed that the extent to which protons are 
polarized is the same in each spin configuration. 
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Figure 7.5 Number of e~^e~ pairs with invariant masses 
Af€[2.7, 3.4] GeV/c^ divided by number of counts in the 
BBC integrated over individual stores, plotted against store 
number, and fit with a constant. 



Eq. 7.2 works in situations where there are enough statistics in A^^t^ N''^, N'^', and N'^'^ 
to determine an asymmetry for a time-period in which Ri, R2, and -R3 are stable. In the 
muon analyses this is possible, and each R is determined for every (approximately eight hour 
long) store. An asymmetry is then calculated by fitting a constant to the asymmetries from 
all stores. In the central analysis there are significantly fewer J/ip counts; so we first assume 



that double spin asymmetries are small so that Eq. 7.2 can be written as 

f N^ - UN^ 



A 



N 



(7.3) 



where the spin orientation of the polarized beam is denoted by the arrow, and there is only 
one relative luminosity IZ = yi- To further simplify, we can eliminate the explicit use of the 
relative luminosity by including both the asymmetry on the left and righljj as 



A 



N 



V 



yivx+yiv^ 



(7.4) 



The geometric scale factor 



/' 



/J'g((/>)sin(^d(/) j^ 
Jo ^ 



R 



I sm 



/f ^('^)rf<^ 



(7.5) 



^Recall from Section Bl that Eq. |4.l| and hence all Ajv formulae above (including Eq. 7.3 1 apply only to the 
left side of the polarized proton. For the asymmetry on the right, an overall negative sign is required. 
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Figure 7.6 On the top, transverse momentum (left) and Feynman-x (right), 
and on the bottom rapidity of accepted e+e" pairs. 



is shghtly different from the factor / in Eq. 7.3 (defined in Eq. 4.2), because left and right are 



treated simultaneously, but this only leads to differences between Eq. |7.3| and Eq. 7.4 of order 



{Ar,f. 



Unfortunately, J/tp statistics in the central spectrometer are still limited to the extent that 



determining asymmetries with Eq. 7.4 from each store is not possible. Instead, statistics are 
integrated across all stores, and a single A^ is calculated for the entire running period. The 
integration across n time periods means that the asymmetry on the left is actually 



^N — - } ^ ^n 1 r ~ ^ / ^^N,i- 



(7.6) 



i=l 



(using Eq. 7.3 for simplicity). Since the measurement is statistically limited, we do not use 
individual fi, Vi, and TZi but instead assume that 



A 



N 



n 

if) Y^ [nJ - in) Nf) 

i=l 

n 



(7.7) 
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is equivalent to Eq. 7.6, where the brackets denote luminosity weighted averages over the course 
of the running period. In order for these two expressions to be truly equivalent it must be true 
that Ui = (7^), fi = if), and Vi = {V) for aU i. 

The acceptance of the central spectrometer is stable enough in the 2006 data that the 
assumption is valid for /. Likewise, variations in the polarization are consistent with systematic 
errors on the measurement of the polarization, meaning that our inclusion of polarization 
uncertainties will take care of any difference between Vi and {V) . Unlike variations in / and 
V, which are under control and can only affect the overall scale of the asymmetry, variations 
in TZ are significant over the course of the measurement and can lead to false asymmetries. 

In order to stabilize 7^ for the measurement in the central spectrometer, a procedure was 
developed wherein several bunches of colliding protons are removed from the analysis so that 
each TZi (the relative luminosity in a given store) is brought as close to unity as possible. 
In principle, any constant value for the relative luminosity will work, but unity is chosen for 
convenience. 

First, all bunches with luminosities greater than 2a away from the mean bunch luminosity 
of a store are removed. A bunch is then chosen at random, and if removing that bunch from 
the analysis brings the relative luminosity closer to unity, it is removed. Otherwise it is kept, 
and another bunch is chosen at random. The process continues until the relative luminosity 
is within 1% of unity or as close to unity as possible given the finite number of bunches. The 
corrected relative luminosities are distributed with an RMS of approximately 1.5% away from 
unity, and the entire procedure removed approximately 5% of the provided luminosity from 
the data sample used in the analysis. 

Since removed bunches are chosen at random, the result of the analysis is not unique, and 
a systematic uncertainty is introduced. In order to determine this systematic uncertainty, 
the analysis is run 5k times, and the resulting asymmetries and statistical uncertainties are 
histogrammed. The mean values of these two histograms are taken as the central value and 
statistical uncertainty in the data point, and the RMS of the histogrammed asymmetry is 
taken as a source of systematic uncertainty. 
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Statistics in the muon analyses are plentiful enough to allow for the measurement of a 
background asymmetry, Afj^ , which makes possible a calculation of the signal asymmetry, 
Aj-j , which is not diluted by background. The signal asymmetry is found with 



A Incl „ A BG 

T 



AJi^ = "^^ ^ ^;"^^ , (7.8) 



where r is the background fraction (Table 7.1). The background asymmetries are determined 



for the analysis of the 2006 dataset using oppositely-charged muon pairs in the invariant mass 
range 1.8 < m{GeY/c^) < 2.5 along with charged pairs of the same sign in invariant mass range 



1.8 < m(GeV/c^) < 3.6 (shaded areas in Fig. 7.3 ^a) and (b) respectively). For the analysis of 
the 2008 dataset the lower limit of the mass range is 2.0 GeV/c^. Background asymmetries 



from this study are shown in Table 7.2 In the central spectrometer, statistics do not allow for 



the measurement of a background asymmetry, and AfP is assumed to be zero in Eq. 7.8 



Table 7.2 Background asymmetries as a function of pT for the PHENIX 
muon spectrometers. 



PT (GeV/c) < xp > dataset Af/^ 

(M3 -0.07 2006 -0.003±0.028 

-0.07 2008 -0.072±0.034 

0.08 2006 -0.008±0.028 

0.08 2008 -0.003±0.035 

0-1.4 -0.07 2006 -0.002±0.031 

-0.07 2008 -0.043±0.039 

0.08 2006 -0.021±0.038 

0.08 2008 -0.060ifc0.046 

1.4-6 -0.08 2006 -0.066±0.050 

-0.08 2008 0.047±0.064 

0.08 2006 0.039±0.056 

0.08 2008 -0.072±0.070 



A final additional complication is that even and odd numbered crossings are triggered by 
separate circuits in the ERT with slightly different gains. This means that detector efficiencies 
can vary depending on the spin orientation and potentially lead to false asymmetries in the 
central spectrometer. We restore the equality in efficiency by measuring asymmetries separately 
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for even and odd crossings and combining the resulting asymmetries. 

7.4 Polarization and Geometric Scale Factors 

Average beam polarizations are measured using the RHIC polarimeters discussed in Sec- 
tion [6ll] and found for the 2006 run to be 



0.53 ± 0.02^^^*^* (clockwise) 

0.52 ± 0.02"^*^* (counterclockwise) 

with standard deviations of 0.03 and 0.04 respectively. The labels represent the direction of 
the beam circulation, and the systematic uncertainties are uncorrelated between beams. There 
is an additional systematic uncertainty of 3.5% correlated between the two beams. 
During the 2008 run, the average beam polarizations were 

0.48 ± 0.02"^*^* (clockwise) 

0.41 ± 0.02''y''* (counterclockwise) 

each with a standard deviation of 0.04 and an additional systematic uncertainty of 3.9% cor- 
related between the beams. 

The deviation of the polarization from nominal at the PHENIX interaction region is mea- 



sured using the local polarimeter discussed in Section 6.2.2 and is found, in radians, to be 



0.064 ± 0.040*^*^* ± 0.086'^y''* (clockwise) 

0.109 ± 0.038*^*^* ± 0.036'^y''* (counterclockwise) 



for the 2006 data and 



0.263 ± 0.030'*''* ± 0.090'y"* (clockwise) 
0.019 ± 0.048*^*^* ± 0.103'^y''* (couterclockwise) 

for the 2008 data. The polarization directions of the counterclockwise-going beam in 2006 and 
clockwise-going beam in 2008 are considered significant and used in the calculation of Aj^. 



74 



The geometric scale factor from 2006 data in the muon spectrometers of / = 1.57 it 0.04 
is determined from J /ij: azimuthal distributions in data and is found to be independent of pT 
within statistical uncertainties. For the 2008 data, the factors are / = 1.64 it 0.01 for the 
clockwise circulating beam and / = 1.56 it 0.01 for the counter-clockwise circulating beam. 
In the central spectrometer, the factor /' is determined from a GEANT (67) Monte Carlo 
simulation of single J/tp decays with a full geometric description of the detector including all 



known inefficiencies, and the resulting geometric scale factors are listed in Table 7.3 



Table 7.3 Geometric scale factors for the central spectrometer determined 
with a Monte Carlo simulation of the detector. 



PT (GeV/c) 


/' 


0-6 


1.62±0.01 


0-1.4 


i.eiibo.oi 


1.4-6 


1.70±0.02 



7.5 Systematic Uncertainties 

Systematic uncertainties can generally be classified as one of three types: uncertainties 
which are uncorrelated between data points (Type A), uncertainties which are correlated be- 
tween data points, usually with an unknown correlation matrix (Type B), and uncertainties in 
the vertical scale which move all data points in the same direction and by the same magnitude 
while maintaining the statistical significance from zero (Type C). 



As discussed in Section 7.3, a Type A systematic uncertainty exists in the central mea- 
surement due to the stabilization of the relative luminosity. There is also Type A systematic 
uncertainty in the muon measurement due to the fit which determines the background fraction 



r in Eq. 7.8 , This uncertainty is estimated by calculating r using a slightly wider mass window 



in the fit to Eq. 7.1 The difference between the resulting A^ and the nominal value is taken 
as the magnitude of the uncertainty. 

Type B systematic uncertainties arise from uncertainties on the polarization which are ei- 
ther uncorrelated between beams or uncorrelated between year and are all 0.003 or smaller 
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(numerical values can be found in Appendix [C]) |j Unlike typical Type B systematic uncer- 
tainties, the uncertainties in this analysis maintain the significance of the data points from 
zero. 

Type C systematic uncertainties come from sources of polarization uncertainty which are 
fully correlated between beams and years (e.g. systematic uncertainty in the hydrogen jet 
measurement due to contamination of the jet). This leads to a scale uncertainty of 3.4% for 
2006, 3.0% for 2008, and 2.4% for the combined 2006 and 2008 datasets. 



A summary of all uncertainties can be found in Table C.l of Appendix^ 



7.6 Results and Discussion 

The final asymmetries included 539ib25 J/tjj mesons in the central spectrometer, 3507ib59 
(3354±58) in the North (South) muon spectrometer from the 2006 run, and 8540±92 (8105±90) 
in the North (South) muon spectrometer from the 2008 run. The values of An obtained by 



analyses from both the central and muon spectrometers are shown in Fig. 7.7 and 7.8, and 



the data points are given in Table C.l of Appendix [Q It is quite exciting to note that the 
Pt integrated data point at (a; i?) =0.08 is 3.3cr less than zero, meaning that there is ~99.9% 
probability that A^ is negative at that xp- The central and backward data points in xp are 
consistent with zero, and there is no observable pT dependence. 

Because there appears to be a non-zero asymmetry at forward xp^ it is interesting to com- 
pare this result with other analyses from PHENIX. Preliminary measurements of the transverse 



SSA for mid-rapidity vr production from PHENIX shown back in Fig. 4.5 cover an x range 
of approximately [0.02,0.3] and show no sign of a non-zero Sivers or multi-parton correlation 
function. The x range covered for the parton in the transversely polarized proton for the xf>Q 



data point in this analysis is shown along with the x coverage for tt production as a function 



of Pt in Fig. 7.9 For forward J /%[) production, the PHENIX muon arms cover approximately 



xi g[0.05,0.15]. 



''These uncertainties would naiVely be Type C, but the different statistical weights used to average the 
uncertainty between beams or years for each data point mean that the magnitude of the uncertainty is not 
necessarily uniform. 
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While the overlap in coverage for forward J/ijj production and central vr^ production at 
Pt ^5 GeV/c mean that the two results appear to disagree, the processes contributing to the 



production of the two particles are also important. Looking back at Fig. 4.6 we can see that at 
Pt ^5 GeV/c, vr*^ production at central rapidity is dominated not by gluon-gluon scattering but 
instead by quark-gluon scattering. This means that any effect from a non-zero gluon Sivers or 
trigluon correlation function for mid-rapidity vr*^ production will be diluted by collisions where 
a quark, not a gluon, enters the scattering from the polarized proton. For J/ip production, 
on the other hand, gluon-gluon scattering dominates for the entire pT range so that there is 
no dilution. Likewise, color factors from the contributing diagrams have been shown to have 
a large effect on the resulting transverse SSA (48) so that a full comparison between the J /tp 
and TT^ transverse SSA will require further theoretical activity. 

At the current state of theoretical understanding it is unclear whether the sign of Aj^i 
is be determined by the Sivers effect or the production mechanism. It is also unclear how 
the magnitude compares to that expected by theory, as (151) gives a possible source for the 
asymmetry but no numerical guidanceFI However, the fact that the J /ijj transverse SSA is non- 
zero in any kinematic region is quite exciting because of the implication that either the Sivers 
or trigluon correlation functions are non-zero. In the TMD approach, this would additionally 
imply that color-singlet diagrams contribute substantially to the J /ip cross section for p + p 
colhsions at -y/s = 200 GeV. 

A great deal of theoretical development and further measurements are necessary to untangle 
both the gluon Sivers effect and the J/ip production mechanism in transverse SSAs (this will 
be discussed further in Chapter^. Nevertheless, the non-zero effect at forward xp which is 



implied by Fig. 7.8 is an exciting and important step. 



^Admittedly, numerical guidance is difficult to obtain, because neither the size of the gluon Sivers function 
nor the J/tp production mechanism have been determined. 
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Figure 7.7 J/ip transverse SSA measured at PHENIX from ^/s =200GeV 
p+p collisions plotted against transverse momentum. Type A 
systematic uncertainties are added in quadrature. Type B sys- 
tematic uncertainties are all less than 0.003 and are not shown. 
Type C systematic uncertainties are also not shown but are 
3.4% (3.0%) for the 2006 (2008) dataset, and 2.4% for the 
2006+2008 dataset. 
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Figure 7.8 J/tp transverse SSA measured at PHENIX from ^/s =200GeV 
p+p collisions plotted against Feynman-x. Type A system- 
atic uncertainties are added in quadrature. Type B system- 
atic uncertainties are all less than 0.003 and are not shown. 
Type C systematic uncertainties are also not shown but are 
3.4% (3.0%) for the 2006 (2008) dataset, and 2.4% for the 
2006-1-2008 dataset. On the left, xp > data points from the 
2006 and 2008 runs are shown separately, and on the right they 
are combined. 
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x1CI= 



^ x10 




Figure 7.9 x range covered for various px ranges of vr'^ production 
at mid-rapidity in the PHENIX central spectrometers (top, 
from (12)) and for the parton in the polarized proton from J/ijj 
production in the PHENIX muon spectrometers at xf>^ from 
p+p collisions at -v/s=200 GeV (bottom). 
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CHAPTER 8. J/i) Angular Decay Coefficients 



8.1 Introduction 

In 2009 the PHENIX experiment collected data from p+p collisions at -^=500 GeV, pro- 
viding access for the first time to J /tp angular decay coefficients from p+p collisions at that 
energy. In this chapter, I will discuss a measurement of the angular decay coefficients A^, X^p, 
and A^^ using the PHENIX muon spectrometers. 

8.1.1 Measuring A^, A,^, and A^^ 

If we require the z-axis of our reference frame to be in the production plane, the full angular 



distribution of leptons from J/f/^ decays is given by Eq. 5.11 It is common, however, to use 



distributions integrated of (p and cos •& written, respectively, as 

ocl + A^cos^i? (8.1) 



dN . . , 2 



d{cos'd) 

and 

dN 2A<o , , 

— «l + --^cos2(^ 8.2 

dip o + A^ 

as expressions for determining A^ and A;^. To determine A^<^ one can then define an angle ip^, 
as proposed in (86), with 

(/? — f for cos t? > 



'P'd 



.3) 



(/? - ^ for cos i9 < 



in order to get a distribution 



dN V2X^ 

-— ex 1 + —— ^ cos 99^ (8.4) 

dip^ 6 + M 



for A^ip. In order to determine the angular coefficients, each of these distributions is measured 
in real data and a simulation with an isotropic decay distribution (A^=A^=A^^=0). To correct 
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for the acceptance, the distributions from data are divided by those from simulation. These 
acceptance-corrected distributions are then fit with Eq. 8.1[ 8.2 and 8.4 to determine each of 
the coefficients. 

One must be very careful, however, to estimate the effect of acceptance on these integrated 
distributions. If the acceptance creates 'holes' in the cos-d-ip phase-space, the acceptance of 
these 'hole' regions cannot be corrected through division by a simulated distribution. The 
relative acceptance of the PHENIX muon spectrometers can be found in Appendix |B] and 
indeed creates large holes. To avoid the effects of holes in the muon acceptance and to include 
the effects of correlations between the various coefficients, we will instead perform a two- 



dimensional fit with Eq. 5.11 To correct for acceptance, we perform a GEANT (67) Monte 
Carlo simulation using a full description of the PHENIX detector to determine the shape of 
the angular distributions in our detectors with an isotropic decay distribution. The angular 



distributions from real data are fit with a convolution of Eq. |5.11 and the simulated acceptance. 



Avoiding muons with momenta near the trigger threshold is quite important, even if the 
simulation exactly describes the trigger. Because the angular distribution is very different for 
A^=l than for A^=0 around the trigger threshold, the physical acceptance of the detector in 
cos 'd can be different for these two situations such that dividing the data by a simulation with 
A^=0 does not properly correct the data for acceptance. An extreme example would be as 
follows: 

Consider two mesons with equal momenta. For the first meson, decay muons are produced 
along the meson's momentum, while, for the second, muons are produced transverse to the 
meson's momentum. One decay muon from the first meson will get a boost along the meson's 
momentum direction, while the other gets a boost in the opposite direction. For the second 
meson, both muons will get a boost in a direction transverse to the meson's momentum. It 
is less likely that both muons from the first meson will make it through the //ID, because 
one of them is being boosted in the wrong direction. Because we require that both muons 
make it through the detector, the trigger turn-on curve will be shifted to higher p^ for the 
first meson. Our acceptance correction will not account for the difference because it uses an 
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isotropic simulated distribution. 

The effect can be seen for realistic angular distributions from A^=0 and A^=l for the 



Helicity frame in Fig. 8.1 where the trigger turn-on curve is clearly shifted to larger pz for 
A^=l. The cut placed on the trigger threshold, pz measured at the first /xTr station, has been 
made tight enough in this analysis to exclude any data around the trigger turn-on. 
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Figure 8.1 pz at the first ^Tr station from a full simulation with A^=l in 
the Helicity frame divided by a full simulation with A^=0 for 
the North (left) and South (right). The A^=l distribution is 
clearly shifted to larger pz Dashed lines show the cut placed on 
this quantity in the analysis. 



8.1.2 Simulation Procedure 

Because we are using a simulation to obtain the acceptance, the simulation itself needs to 
be well tuned to reproduce both low-level cluster-related quantities and high level kinematic 
distributions. In this section, I will give a general outline of the procedure used to tune the 
simulation to data. More detailed information and plots from the simulation will be given 
throughout the rest of this document. The procedure will be as follows (a fiowchart can also 



be found in Fig. 8.2): 



1. Introduce inactive areas into the simulation. 

2. Tune the total cluster charge, tracking resolution, and cluster shape to match those found 
in real data. 
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3. Use the tuned simulation to determine the acceptance x efficiency and yield for J/tp 
mesons as a function of transverse momentum {px) and rapidity (y). 

4. Generate a second set of events using the shape of the measured yield in pT and y from 
Step 3 and an isotropic decay distribution. Use this simulation to determine A^, A,^, and 

5. Repeat Step 3 using the measured A^, X^p, and A^;^ from Step 4. Continue iterating until 
the shape used in Step 3 is consistent with the shape of the measured yield. 



Remove inactive 

detector areas from the 

simulation 



T 



Tune simulated clusters 
to match data 



r 



Calculate 
acceptance x efficiency 

and yield using 
measured X^,X , and X^ 



X 



Generate simulated 

angular distributions 

with X =X =X =0 

8 ,p 3,p 



X 



Detemiine X^,X , and X 



Iterate until measured 
shape in p and y is 

consistent with the 
simulated shape 



Figure 8.2 Flowchart of the iterative procedure for tuning the simulation 
to real data. 



8.2 Quality Assurance 



Before the data can be taken seriously, we need to place criteria on both the subset of the 
data used in the analysis and the tracks in the detector considered to be muons. A deter- 
mination of the angular decay coefficients is especially sensitive to large changes in detector 
acceptance; so we will use criteria which require that the acceptance stay relatively constant 
over the course of the analysis. Requirements will also be placed on tracks to increase the 
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number of muons from the J/ip signal relative to background from random track combinations 
in the detector. 

8.2.1 Run Selection 

Data is taken by the PHENIX experiment in approximately one hour segments (called 
runs) to allow for the rejection of data which are taken with less than ideal conditions. To 
ensure that the acceptance is approximately constant across all runs in the analysis, runs are 
required to meet the following criteria: 



< 4(59) tripped or disabled high voltage channels in the North (South) /xTr (see Fig. 8.3). 
Each channel corresponds to approximately 1/4 of one octant in a gap at one of the three 
stationso 



• < 2(3) tripped high voltage channels in the North (South) /xID (see Fig. 8.4). Each 
channel corresponds to approximately 1/3 of one of the 6 panels which make up either 
the horizontal or vertical plane of one gap. 

• < 1% of the data packets lost from the /iTr during transmission to the data acquisition 
system. Losing one packet of data corresponds to a loss of that event for approximately 
one octant at one station in the /xTrjj 

• runs should not have been ended due to problems with any of the muon detectors 

• the data acquisition system should indicate that the run ended cleanly 

Of the 259 runs applicable to the analysis, the run selection criteria left 236 (238) runs for anal- 
ysis in the North (South) muon spectrometer. The total integrated luminosity was 17.6 pb^^, 
and after run selection we are left with 16.0 (16.5) pb~^ in the North (South) spectrometer. 



'^There are significantly more disabled channels in the South ^Tr than the North because of high backgrounds 
in the inner channels of the first station. 

^There are only two data packets for the each pID; so we require that all /iID data packets be received. 
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Figure 8.3 Number of tripped or disabled channels in the /iTr his- 
togrammed by run for the arm on the North (South) side of 
PHENIX on the top (bottom). 
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Figure 8.4 Number of tripped or disabled channels in the /xID his- 
togrammed by run for the arm on the North (South) side of 
PHENIX on the top (bottom). 



8.2.2 J/tp Selection 

Criteria for selecting J/ip mesons must simultaneously reduce background and maintain 
signal. In order to determine whether or not a pair of tracks come from a J/il^, we need to 
ensure that both tracks behave like muons and originate from a real collision. 

For an event to be considered, the following conditions must be met in either the North or 
South muon spectrometer: 

/xID 2-deep and no vertex BBC Trigger Two deep requirements must be met in the ei- 



ther the North or South /ilD LLl, discussed in Section 6.2.5, along with at least one tube 
fired in both the North and South BBC. 



BBC z-vertex The position of the event vertex along the beam direction must be within 



40 cm of the center of the central spectrometer as measured by the BBC in order to 
avoid steel from the magnets located at ~45 cm. 

DDGO. The angle between the track in the /iTr and the track in the //Id projected to the 
closest /ild gap to the interaction point must be less than 9° to ensure that tracks in the 
two detectors are associated with each other. 

Vertex y^ /ndf A fit of the two tracks in the /xTr to the vertex measured by the BBC must 
have a x^ /ndf <4 for 4 degrees of freedom. Such a requirement removes accidental 
combinations of random tracks when constructing the pair. 

Minimum pz at the first //Tr station The magnitude of the track momentum along the 
beam direction measured at the first /iTr Station must be greater than 1.75 (1-45) GeV/c 
in the North (South) spectrometer to avoid the turn-on curve associated with the trigger 



(Fig. 8.1). 



8.3 Low Level Simulation Tuning 

To ensure that our simulation correctly reproduces data, we need to first check basic quan- 
tities like the size and shape of charge clusters in the yuTr as well as dead areas in both the 
/iTr and /ilD. 

Charge is deposited in a cluster of strips by particles traversing the //Tr, and the total 
charge in those clusters is simulated by a Landau distribution 

„ /■CXD 

P{q) = go + - / e-*'°s*-^*sm(7rt)(it (8.5) 

where the parameters go and gi are chosen octant-by-octant for each gap in the detector such 
that the distribution of total cluster charge matches the data. The total cluster charge from a 



simulation of single muons is compared with cluster charges from data in Fig. 8.5 and 8.6 and 
agrees quite well after the tuning of go and gi . 

The simulated shape of clusters in the /xTr is determined by a Mathieson distribution (91; 
131). The capacitive coupling between neighboring cathodes and the spacing between cathode 
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Figure 8.5 Comparison of total cluster charge in the North //Tr in simula- 
tion (open circles) and data (closed circles). 
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Figure 8.6 Comparison of total cluster charge in the South /xTr in simula- 
tion (open circles) and data (closed circles). 
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Figure 8.7 Comparison of peak strip charge over total cluster charge in the 
North yuTr in simulation (open circles) and data (closed circles). 



and anode are tuned for each octant in each gap so that the fraction of the total cluster charge 
carried by the peak strip matches between simulation and data, and a comparison between the 



simulated and measured distributions is shown in Fig. 8.7 and 8.8, While these distributions do 



not match perfectly, the effect of the cluster shape on acceptance should be minimal compared 
to differences in the total charge, which couple directly to the efficiency of the detector. 

Previous PHENIX analyses have assumed that both the total cluster charge and cluster 
shape are constant across a given gap in the /iTr. Through the course of this analysis, I found 
that there are significant variations from octant-to-octant within a gap (due to differences in 
the electronics which measure gain). The analysis presented here is the first to use simulated 
cluster charges and shapes which are not only tuned gap- by-gap but also octant-by-octant. 

In addition to the size and shape of clusters, it is important that the tracking resolution 
of the simulated ^Tr match the resolution of the data. This resolution is measured by the 
'w' coordinate, defined as the distance along a direction perpendicular to the strips in the 
chamber. Unfortunately, the resolution of the simulated ^Tr is typically much better than the 
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Figure 8.8 Comparison of peak strip charge over total cluster charge in the 
South /uTr in simulation (open circles) and data (closed circles). 



real detectorjj and in order to correct for this difference, the simulated charge distributions 
are smeared by an additional amount A 



Qsmeared 
^Qsmeared 



Qgenerated i Ji-T 



A- 



9 



(8.6) 
(8.7) 



where P is a random number chosen from a normal distribution centered about the origin with 
<T=1, r and g are the root mean square ADC counts (rms) and gain, respectively, measured 
from the calibration of the strip with charge q. The parameters A for each /xTr gap are tuned 
so that the simulated difference between the w coordinate of the cluster and the track (the 
residual) match the measured distributions. The simulated and measured tracking residuals 



are shown in Fig. 8.9 and 8.10 



In addition to the cluster charge distributions and residuals, it is important that the physical 
acceptance of both the /iTr and fjJD match between simulation and data. For the fjJD this 
matching is accomplished by measuring the efficiency of each set of two larocci tubes in both 



^The diflference is likely due to misalignments in the real detector which are not properly accounted for in 
the reconstruction or simulation. For example, the large width of the residual in the North Arm Station 0, Gap 



seen in Fig. 8.9 is entirely due to a 0.1 cm misalignment found in that gap after the data was produced. 
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Figure 8.9 Comparison of the tracking residual A-u;, the difference between 
the cluster and track w coordinate in the North fxTr between 
simulation (open circles) and data (closed circles). 
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Figure 8.10 Comparison of the tracking residual Ait;, the difference be- 
tween the cluster and track w coordinate in the North /iTr 
between simulation (open circles) and data (closed circles). 
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the North and South detectors and applying those efficiencies to the simulation. The //ID 
efficiency was determined from data using tracks with hits in nearly every 2-pack (pair of 
tubes). The efficiency is calculated as 

Tracks with hits in all 10 planes 
p ano o m eres Xracks with hits in at least 9 planes, excluding the plane of interest 

and is projected to 2-packs in the /ilD to determine an efficiency for each 2-pack. The average 
efficiency is found to decrease with interaction rate (Fig. 8.1l[ )p]but we simply use the mean 
value of the efficiency over the course of the running period for each 2-pack, as a uniform 
change in efficiency should have little effect on the acceptance to angular distributions. The 



mean is calculated by integrating both the numerator and denominator of Eq. 8.8 over the 
course of the running period. 



I'lntif 1*66/27 

Constant .91 85 ± 0.0005251 
Slope -D.002076±4.104b-05 




^ q_^4^1aXaaaXaaaXaaaXaaa^aaaXaaaXaaaXaaaXaaaXaaaXa 



o 0.93t- 

m 



X^/ndf 407.2/27 

Constant 0.9443 ± 0.0005063 
Slope -0.0032351 3.81 1e-0S 
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Figure 8.11 Average efficiency for pairs of larocci tubes (called '2-packs') 
in the //ID plotted against BBC rate in the North (South) arm 
on the left (right). 



In the //Tr, channels are removed from the simulation if the high voltage for that channel 
was off for the majority of the running period. Fiducial cuts are then made to the dead regions 
in order to ensure that no tracks reconstructing into those regions are used. 

After low-level tuning and fiducial cuts are applied, we can look at the distributions of the 
azimuthal and radial track projections on each /iTr station to see that the simulation adequately 



''The decrease in efBciency with interaction rate is due mostly to an associated increase in beam-related 
backgrounds. As the number of hits in the /ilD increases, the average current in each tube is increased, 
decreasing the voltage across the tube (which decreases the efficiency). 
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reproduces the data. These distributions are shown in Fig. 8.12 and 8.13 respectively and show 
good agreement between the simulated distributions and those from data. 
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Figure 8.12 Comparison of azimuthal track projections between simulation 
(open circles) and data (closed circles) in the North (left) and 
South (right). 



8.4 J/ip Yield 

A necessary step in measuring the angular decay coefficients for J/^ mesons is to determine 
the shape of the J/ip yield in px and rapidity in order to have a proper kinematic shape for 
the simulation. Since no previous measurements have been made of the J/V' yield from p+p 
collisions at -y/s=500 GeV in 1.2< \y\ <2.2, we will need to determine this shape and the 
uncertainty associated with it from data. 

In our data, the kinematic distributions of /U+/i^ pairs in the J/^ mass region contain some 
background, and it is important to remove the contributions from this background. Two types 
of background contribute to the number of oppositely charged muon pairs in the J/ip mass re- 
gion: (1) random combinations of tracks in the detector (called the combinatorial background), 
and (2) physical processes like p+p^T>DX^fj,^ fi^X (open charm), p+p^JiBX^fi^fi^X 
(open bottom), and Drell-Yan production, which all produce a continuum of correlated ^^ ^~ 
pairs (called the continuum background). 

To estimate the combinatorial background we will us the same method discussed in Sec- 
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Figure 8.13 Comparison of radial track projections between simulation 
(open circles) and data (closed circles) for each /^Tr station 
in the North (left) and South (right). 



tion 



7.2, which is to take twice the geometric mean between the number of positive and nega- 



tively charged muon pairs, 2^/N^-i 



+ 7V - - 



. Physical background is more difficult to estimate 
than combinatorial background, but we are helped by the fact that the physical backgrounds 
form a continuum which can be approximated with a falling exponential under the J/t/j peakr] 
The following functional form is used to describe the dimuon continuum in mass: 



dN 
dM 



Ae 



-bM 



N 



+ 



J/V 



("-^^J/,/.)^ . ,. (M-Mj^^)'^ 



2^^J/^,i 



J/V),2 _|_ 



f2Nj/^ 



J/i',2 



(8.9) 



The second normal distribution accounts for large tails in the J/'0 mass peak found in simula- 



""While there is no clear theoretical reason that an exponential should describe the combination of the Drell- 
Yan and open-heavy flavor backgrounds in the continuum, such a shape has historically parameterized the 
background very well. 
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tion, and the /2 is fixed by simulation. The J /ip resonance is quite wide in data because of the 
resolution of the spectrometer. Because of this width, the ip' resonance is not very prominent, 
especially in regions of little acceptance, and is not included in the fit. 

The exponential shape of the background is modified at low mass by our trigger require- 
ment. In order to take this into account and obtain a better fit, a simple Monte Carlo simulation 
is used to determine the effects of the acceptance and trigger on an exponentially falling back- 
ground. Lepton pairs are generated with a normal distribution in rapidity and an exponentially 
falling distribution in mass and are put through a simple model of the detector acceptance and 
trigger requirement. The output distributions from simulation are then applied as a weight to 
the exponential used in the fit in order to obtain a more realistic distributionjj 
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Figure 8.14 Example plots of the total counts of oppositely charged muons 
from one kinematic bin along with the associated combinato- 
rial background estimation (left). Example fits to the back- 
ground subtracted distribution (right). Separate fit ranges are 
show in red and black, and dashed lines show the mass range 
in used to determine the signal. 



Varying the range for the fit to Eq. 8.9 can potentially change the number of J /il) mea- 
sured. In order to account for this variation, fits are repeated with two separate fit ranges: 
Mg[1. 8,7.0] GeV/c^ and Mg[2. 2,6.0] GeV/c^. Separate fits are also done with /2 increased and 
decreased by 25%. 



''Another way of taking t he e ffect of the trigger into account is to use a polynomial to describe the background 
shape, as we did in Section 7.2 
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The separate fit ranges and 3 scenarios for /2 account for 6 fits, and they are plotted 



together in Fig. 8.14, where the 3 separate scenarios for /2 are indistinguishable. The mean 
of the integrals of the normal J/V' distributions from the 6 signal extractions is taken as 
the central value of the data point for each bin in pT and rapidity, and the RMS about the 
mean number of J/V' from all signal extractions is taken as a source of systematic uncertainty 
which is mostly uncorrelated between bins. We will assume that this uncertainty is entirely 
uncorrelated between bins and add it in quadrature with the statistical uncertainty. 

When fitting the yields in pT and rapidity, it is important that the data points be located 
at the mean px or rapidity for that bin (depending on which is being histogrammed). To 
determine the mean px and rapidity for each bin, we generate a sample of simulated J /tp 
mesons and put them through a GEANT (67) model of the detector, keeping track of both 
the generated and reconstructed kinematics. For each bin in reconstructed pT and rapidity we 
histogram the generated quantities, taking the mean value as the centroid of the bin in our 
final yield and the RMS as the uncertainty for that bin in abscissa. 

After the data points have been determined, we need to parameterize the J /ij) kinematics 
for use in simulation. A Kaplan function (119; 114) is used to parameterize the px distribution 

' '"-4l+(?^VV", (8.10) 



iTTprdpT y \PoJ ^ 

and a normal distribution centered about zero parameterizes the shape in rapidity. To get a 
reasonable distribution of event vertices in our simulation, vertices are taken from a sample of 
events for which at least one tube fired in both the North and South BBC. 



As discussed in Section 8.1.1, the method for determining the J/^ yield is iterative. First 
we measure the yield using the acceptance x efficiency from a given simulation, then we 
repeat using the measured yield shape as input to the simulation. The final iteration used 
to determine the acceptance x efficiency had n=4.5 and po=3.4: in the Kaplan function and 
a=1.27 for the normal distribution in rapidity. For the angular decay coefficients we use a 
linear interpolation between the measured angular decay coefficients of the Helicity frame as a 



function oi px (consistent with the final measured values). Fig. 8.15 shows the acceptance x 



efficiency from the simulation plotted against px and rapidity, and Fig. 8.16 shows the yields. 
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From the fit parameters shown in the latter figure, it is clear that the input shape in pT and 
rapidity is consistent with the measured shape in both spectrometers. 

Because the parameterizations of the yield shape in px and rapidity are not perfectly 
determined, we will need to introduce systematic uncertainties to account for the full range of 
fits in both the North and South muon spectrometers. When determining the uncertainty, a 
10% systematic uncertainty, correlated between data points, is attributed to the acceptance x 
efficiency (because our simulation is not an exact description of the detector) . Yields are fit with 
a modified likelihood (described in (11)) which takes these correlated systematic uncertainties 
into account. The resulting parameters from the fit are then averaged, and the values at itlo" 
are taken as the extreme shapes of the yield, used to determine the systematic uncertainty in 



Section 8.5, The nominal Kaplan function in pT was found to have n = 4.52 and po = 3.41, 
while the steep function at the la limit had n = 4.43 and po = 3.46 and the shallow function 
had n = 4.61 and po = 3.35. The nominal normal distribution in rapidity had a = 1.29, while 
the steepest and shallowest distributions had a = 1.25 and a = 1.33 respectively. 
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Figure 8.15 Acceptance x efficiency plotted against transverse momentum 
on the left and rapidity on the right for the North (open circles) 
and South (closed circles) muon spectrometer. 
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Figure 8.16 Acceptance corrected J/tp yield plotted against transverse mo- 
mentum on the left and rapidity on the right for the North 
(bottom) and South (top) muon spectrometer. Dashed lines 
show the average fits to a Kaplan function on the left and nor- 
mal distribution on the right, and solid lines show the fits at 
iblcj. 



8.5 J/ip Angular Decay CoefRcients 

To determine angular decay coefficients, two dimensional histograms of 5 bins in cos'i? by 8 
bins in if are filled from data and an isotropic simulation (plots of the acceptance in all frames 
are shown in Appendix pi). The combinatorial background is then estimated and subtracted, 



and the mass distribution is fit with Eq. 8.9 as discussed in Section 8.4 to extract the number 
of J/tp mesons in each kinematic bin. Finally, two dimensional histograms of the number of 



J/ip in cosi?-(^ from data are fit with the full angular distribution of Eq. 5.11 convoluted with 
the histogram from simulation. 

We quantify the uncertainty due to the input shape of the yield by using the extreme 



shapes in px and rapidity described in Section 8.4, A cos'd-(p surface is drawn for the nominal 



shape in pT and the two extreme shapes in rapidity as well as the nominal shape in rapidity 
and the two extreme shapes in pT- The surface is then varied 15k times with a Gaussian 
probability about the nominal shape with the extreme shapes at itlcr. The resulting coefficients 
are histogrammed to determine the systematic uncertainty separately for the px and rapidity 
shapes. These uncertainties are summed in quadrature and considered correlated between data 
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pointsjj 

The contribution of background to the angular decay coefficient A^ is 

^■&,S — —^ ^-0,3 5 lo-J^-Lj 

where Rs is the fraction of pairs which come from J/^|J (a derivation can be found in Ap- 
pendix IaI). Likewise, for A^ and A^^ we have 

, _ \^ 3_+jVs_ _ . (1 - Rs) 3 + X^^s (R^'2) 

-Rs 3 + A^^fit ' -Rs 3 + A^3 

and 

■, _ A^y^fit 3 + A,?,s , (1 " Rs) 3 + A^^s fs n^i 

Rs 3 + A^^fit ' -Rs 3 + Atf^B 

(Derivations of these expressions are in Appendix [A|) . In practice, the background coefficients 
A^.B, etc. are very difficult to quantify. Even if the combinatoric background is subtracted 
and the decay coefficients of the continuum are assumed to be negligible, the difference in 
Pt and rapidity shape between the continuum background and signal can cause large false 
asymmetries to be measured. To make matters worse, the Drell-Yan background, which makes 
up a small part of the continuumFjis known to have a large X^ in the Collins-Soper frame (152). 
To remove background contributions, fits are performed in each two-dimensional bin of the 



COS1? — ip distribution using the distribution in Eq. |8.9| as described in Section 8.4 

Measured angular decay coefficients included approximately 13167 J/ip mesons in the North 
and 25390 in the South muon spectrometer. After all fits are performed and systematics 
included for the input kinematic shapes, the resulting angular decay coefficients are plotted 



against px and xp (Fig.|8.17). A single px and xp integrated data point is shown in Fig. 8.18 



and numerical values for each data point can be found in Appendix [C} 

8.6 Discussion 

The integrated A^ coefficient is consistent with zero in all reference frames, but the coeffi- 
cient appears to be negative at low pT in the Gottfried- Jackson forward frame. Unfortunately, 



'^Note that the correlation matrix of this uncertainty between data points is different for angular decay 
coefficients plotted against pr than it is for coefficients plotted against xp. 

*The Drell-Yan contribution is estimated in the J/tp mass region to be approximately 2% of the J/ip yield 
from a PHENfX analysis of p+p colhsions at ^=200 GeV (90) 
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Figure 8.17 From top to bottom: Angular decay coefficients from the 
Collins-Soper, Gottfried-Jackson forward, and Gottfried-Jack- 
son backward, and Helicity frames plotted against pT on the 
left and xp on the right. Point-to-point uncorrelated uncer- 
tainties are shown with error bands, and point-to-point corre- 
lated uncertainties with boxes. 
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predictions do not exist yet for J/^ decay coefficients from p+p collisions at ^/s=500 GeV, but 



the prediction from the CSM at NLO shown back in Fig. 5.9 for p-\-p collisions at ^/s=200 GeV 



appears to be consistent with the measured A^ in the Helicity frame, and the angular decay 
coefficients are not expected to change drastically with ^/s. Predictions from the COM are 
typically valid for a minimum px of 5 GeV/c, larger than the px of this measurement. 

From a phenomenological point of view, the Gottfried-Jackson Forward frame is slightly 
preferred as closest to the 'natural' reference frame for the J/V' spin alignment (the frame 
where only A^ is non-zero) as A^ is largest in magnitude in that frame, while A;^ and A^i^ are 



smallest (Fig. 8.18). Measurements from the HERA-B experiment for p-|-N collisions (3) found 
the Collins-Soper frame to be the closest to natural, but those experiments were carried out 
in a fixed-target environment where only a single Gottfried-Jackson frame is relevant. 
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Figure 8.18 Magnitude of the angular decay coefficients in the Collins- 
Soper (CS), Helicity (HX), Gottfried- Jackson Forward and 
Backward (GJF and GJB) integrated across pT and xp. Statis- 
tical and systematic uncertainties are summed in quadrature. 
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Values of 



It is also interesting to look at the frame-invariant A directly using Eq. 
A along with their uncertainties are plotted for each reference frame against px and xp in 
Also plotted in that figure are the mean values of A assuming that uncertainties are 



Fig. 



8.19 



completely correlated between frames. The measured A are consistent between frames within 
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their systematic uncertainties, and it is quite interesting to note that the mean A is consistent 
with a non-zero longitudinal spin- alignment which increases in magnitude with increasing px- 
This suggest that there exists a natural reference frame wherein the J /ij: has a longitudinal 
spin-alignment and that none of the measured frames are optimal for maximizing | A^ | . 
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Figure 8.19 Frame-invariant angular decay coefficient, A, from measure- 
ments in the Collins-Soper (CS), Helicity (HX), Gottfried- 
Jackson Forward and Backward (GJF and GJB) reference 
frames plotted against px and xp on the top. On the bottom, 
the average of the four measuremnets where wide lines cor- 
respond to uncertainty uncorrelated between points and wide 
boxes correspond to correlated uncertainty. 



To summarize, we can make several strong statements regarding the J /tp angular decay 
coefficients from these measurements: 



A^ is small in each of the four measured reference frames. 



• The Gottfried-Jackson Forward frame appears to be the closest to the natural reference 
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Table 8.1 Frame-invariant angular decay coefficient, A, calculated for 
a single data point integrated over pT G [0, 7] GeV/c and 
\xf\ G [0.015,0.06] 



Reference Frame 


A 


uncorr. 


corr. 


Collins-Soper 


-0.14 


+0.08 
-0.11 


+0.07 
-0.10 


Gottfried- Jackson 
Forward 


-0.25 


+0.06 
-0.08 


+0.02 
-0.04 


Gottfried- Jackson 
Backward 


-0.16 


+0.08 
-0.11 


+0.10 
-0.12 


Hclicity 


-0.17 


+0.09 
-0.11 


+0.10 
-0.13 



frame, with the smallest \\^p\ and |A^<p| and the largest |A^|. 

• Using the frame invariant approach proposed in (86), we find a A which is likely negative 
and increasing in magnitude with increasing pj^. 

All of these conclusions give useful information about the J /ip production mechanism and can 
be used in correlation with future measurements from other collision systems and energies to 
determine how the J /ij) is produced. 
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CHAPTER 9. Conclusions and Future Measurements 



Two measurements have been presented in this document: a transverse single spin asym- 
metry (SSA) of J/'ijj mesons and a comprehensive measurement of all relevant J/V' angular 
decay coefficients in various reference frames. Both measurements show information which is 
new and relevant. 

The transverse SSA of J /ip mesons is 3.3cr less than zero at <Xir>=0.8 and consistent with 
zero for mid and backward xp- The trend of having a non-zero transverse SSA at forward xp is 
consistent with transverse SSA measurements of other particles but is unexpected in the case 
of the J /ip- The fact that the transverse SSA is non-zero potentially implies a non-zero gluon 
Sivers or trigluon correlation function in the proton. It might also mean that forward J /^ip 
production from p+p collisions at -^8=200 GeV is dominated not by color-octet production 
but more likely by color-singlet production. The calculation of trigluon correlation functions 
for J /ijj production would clarify the relationship between the transverse SSA and production 
mechanism. 

Measurements of the J/V' transverse SSA with other center of mass energies and collisions 
systems would help to confirm the existence of a transverse SSA at forward xp and to determine 
the mechanism for such a transverse SSA. While the evolution in Q^ of the gluon Sivers and 
trigluon correlation functions are unknown, transverse SSAs for other particles have persisted 
through a wide range in y/s in various collision systems. A measurement of the J /ip transverse 
SSA in SIDIS would be especially useful, as (151) predicts a vanishing asymmetry for the 
color-singlet model in SIDIS but non-zero asymmetry for the color-octet model. Mapping 
the dependence of the effect on collision system, therefore, could reveal considerably more 
information about the J /ip production mechanism. 
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The J /ip spin-alignment from p+p collisions at y^=500 GeV for 1.2<|y|<2.2 is consistent 
with zero in all measured reference frames when integrated over p^ and xp- The Gottfried- 
Jackson forward frame appears to have a negative A^ at low pT and is the closest to the natural 
reference frame. The frame invariant A implies that there exists a natural frame wherein A^ 
becomes increasingly negative as px increases. 

Measurements of all J /ij: angular decay coefficients (not just A^) using various collisions 
species and energies would lead to a much better understanding of the production mechanism, 
as the kinematic dependencies of the color-singlet and color-octet diagrams are quite different. 
In particular, a measurement of this kind at larger p^ could definitively show whether or not 
the Color-Octet Model describes data, and if not, in what respects it disagrees (i.e. is the 
Helicity frame really the natural frame for J /ip spin-alignment?). 

The J/V' production mechanism provides access to basic QCD dynamics, but the road to 
understanding it has not always been smooth, and we have not yet reached the end. The 
measurements presented in this document represent the application of novel ways of thinking 
about the production mechanism, and I hope that these measurements, along with similar 
measurements from other experiments, will lead to a better understanding of the J /ip in 
particular and QCD in general. 
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APPENDIX A. Background Contributions to Angular Distributions 



This appendix presents a derivation of the background contribution to the angular de- 



cay coefficients discussed in Section 8.5, If we exphcitly include background in Eq. 8.1, the 
differential yield of inclusive pairs Nj becomes 

-^ « /s (1 + A^,scos2^) Aesicos^) + /b (l + A^.bcos^^) Aesicosi^) (A.l) 

where /s and /b are the fraction of pairs in the signal and background respectively and A^^S; 
A^.B Eire the signal and background coefficients. To simplify the expression, we will assume 
that the acceptance times efficiency Ae{cos'&) is the same for the signal and background: 

J^oc(l+ ^'^'f + ^f'^^ cos^ A Aeicos ^) (A.2) 

dcos'd V JS + JB / 

If we define Rs = r fr , we can then see that A^ determined by our fit to the inclusive 
distributions is 

Atf,fit = i?sAtf,s + (1 - i?s)A^,B, (A.3) 

which can be rearranged to give Eq. |8.11 



The derivation for A;^ and A^,^ are very similar, and we will only present the derivation for 



A(p. Explicitly including background in Eq 8.2 we have 



"^^ oc/sfl + -^^cos2(^Ves(^)+/Bfl + ;,^^cos2(^)A6B((/^). (A.4) 



d{(p) V 3 + A^,s / V 3 + A^,B 

Assuming that Aes{if)=AeBi^), this expression simplifies to 

dN ( 2RsX^,s „ ^ 2(1 - Rs)X^,B . , , , , ,, ,, 

oc I 1 + -^^cos2(^H — Y — cos 2ip \ Ae[ip). (A. 5) 



d{ip) V 3 + Atf,s (3 + A^,b) 

from which we can identify 



Xipfit _ RsX<p,s (1 - Rs)\,B / . „x 



3 + A.|y^fit 3 + A^^s (3 + A^^B 
which can be rearranged to give Eq. |8.12 
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APPENDIX B. Relative Acceptance for Decay Angles 
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Figure B.l Relative acceptance for cos-i^-c/? in (from top to bottom) the 
Collins-Soper, Gottfried-Jackson Forward, Gottfried-Jackson 
Backward, and Helicity frames for increasing px from left to 
right. Plots are drawn from real data for the North muon 
spectrometer, and the acceptance is nearly identical for the 
spectrometer in the South. 
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Figure B.2 Relative acceptance for cos'd-ip in (from top to bottom) the 
Collins-Soper, Gottfried-Jackson Forward, Gottfried-Jackson 
Backward, and Helicity frames for increasing xp from left to 
right. Plots are drawn from real data for the North muon 
spectrometer, and the acceptance is nearly identical for the 
spectrometer in the South. 
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APPENDIX C. Data Tables 



This appendix contains a collection of data tables for the J/V' SSA presented in Chapter M 
and the J/^ angular decay coefhcients presented in Chapter l8J 
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